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INTERACTIONS  OF  RADIATION,  MIGROPHYSICS, 

and  turbulence  in  a  two-dimensional 

CIRRUS  CLOUD  MODEL 


Yu  Gu  *  and  K.  N.  Liou 
Department  of  Meteorology/CARSS 
Univeisity  of  Utah,  Salt  Lake  City,  Utah 


ABSTRACT 

A  two-dimensional,  time-dependent  numerical 
model  has  been  developed  for  the  studies  of  the 
interactions  among  turbulence,  radiation,  and  microphysics 
in  cirrus  clouds.  The  model  incorporates  parameterizations 
that  account  for  phase  changes  of  water  and  the  effects  of 
microphvsical  structure  on  the  ice  water  vertical  flux .  Also 
included  is  an  advanced  interactive  radiative  uansfer 
scheme  which  developed  at  University  of  Utah.  To 
investigate  the  interactions  of  various  physical  processes,  a 
complete  version  of  the  second-order  turbulence  closure 
model  has  been  developed.  Numerical  simulation  of  a  thin 
cirrostratus  shows  that  the  radiative  processes  begin  to 
play  an  important  role  once  a  sufficient  amouirt  of  ice  water 
is  produced.  The  induced  turbulence  has  significant  effects 
on  the  phase  changes  and  the  formation  of  ice  \yater. 
Turbulence  also  enhances  the  radiative  net  cooling  in  the 
ice  water  region,  which  is  an  important  mechanism  for  the 
maintenance  of  cilTUS  clouds. 


1.  INTRODUCTION 

Ubiquitous  cirrus  clouds  play  an  important  role  in 
the  radiation  field  of  the  earth-atmosphere  system,  and 
hence  the  earth's  climate.  However,  fundamental 
understanding  of  the  mechanisms  for  their  formation  and 
dissipation  is  still  extremely  limited  and  has  received 
substantial  scientific  attention  of  late  G-iou  1992). 

The  simulation  results  of  Starr  and  Cox  (1985) 
show  that  the  effects  of  radiative  processes  and  vertical 
transports  arc  both  significant  in  the  cirrus  cloud  formation 
and  maintenance.  Other  cirrus  cloud  modeling  studies 
(Jensen  et  al.  1994)  used  simplified  parmeterizations  for 
radiative  processes  and  have  not  considered  turbulence 
effects  in  the  model.  At  this  point,  interactions  amorrg 
radiative  processes,  turbulence  and  ice  microphysics  in 
cirrus  clouds  still  remain  poorly  understood.  In  the  present 
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paper,  we  develop  a  two-dimensional,  time-dependent 
numerical  model  for  the  interaction  studies. 


2.  MODEL  description 


2.1.  Governing  Equations 

The  mean  field  equations  governing  the  modons  in 
the  x-z  plane  based  on  the  Boussinesq  approximation  can 
be  written  as  follows; 

V)  -  Wo  5+ E- <,  -  0+ 


'ai  Oo  _  _ 

3*(u*w’)  3^(u*  w*>  *) 
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(2.2) 


(2.3) 
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where  the  subspripts  refer  to  the  initial  state,  which  is 
horizontally  uniform,  cloud  free  and  at  rert  except  for  a 
constant  vertical  motion  Wq',  the  unsubscripted  yanables 
are  regarded  as  mean  fields  based  on  the  initM  state; 
variables  with  prime  represent  turbulent  fluctuations;  Qr 
and  Qc  are  the  potential  temperamre  tendencies  due  to 
radiative  processes  and  phase  chMges 
respectively;  v  is  the  ice  water  relative  fall  speed;  C  is  me 
subUmation  between  water  vapor  and  ice  water,  q  is  water 
vapor  specific  humidity;  t  is  ice  water  mixing  rauo,  v  is 
eddy  viscosity;  k,,  k,.  k,  are  thermal,  water  vapor  an 
ice  water  eddy  diffusivities,  respectively;  €  is  related  to 
the  ratio  of  molecular  weight  of  water  to  that  of  dry  air,  and 

other  notations  are  conventional. 

The  governing  equations  (2.l)-(2.4)  require  mr^ 
turbulent  fluxes,  which  can  be  predicted  using  the  second- 
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order  turbulence  closure  model.  The  closure  assumptions 
are  described  in  Mellor  and  Yamada  (1974)  and  arc  not 
given  here. 

2.2.  Radiative  Transfer  Parameterizations 

The  radiative  transfer  scheme  we  have  used  is  based 
on  the  dclta-four-stream  approximation  developed  by  Liou 
et  al.  (1988).  It  can  be  applied  to  a  nonhomogeneous 
cloudy  atmosphere  by  dividing  it  into  a  number  of  layers 
each  of  which  can  be  considered  homogeneous  with  respect 
to  the  single-scattering  albedo  and  phase  function.  A 
parameterization  of  non-grtQr  gaseous  absorption  based  on 
the  correlated  k-distribution  method  has  been  developed  by 
Fu  and  Liou  (1992).  In  this  method,  the  cumulative 
probability  g  of  the  absorption  coeffleient  in  a  spectral 
interval  is  used  to  replace  wavenumber  as  an  independent 
variable.  Absorption  due  to  HjO,  O3,  CO2.  and  O2  is 
accounted  for  in  the  solar  spectrum.  In  the  thermal  infrared, 
absorption  due  to  HjO,  Oj,  CO2,  CH4.  and  N2O  is 
included.  Using  the  correlated  k-distribution  method,  121 
spectral  calculations  are  required  for  each  vertical  profile. 

A  new  parameterization  has  been  developed  for  the 
calculations  of  the  single-scattering  properties  for 
hexagonal  ice  crystals.  In  this  parameterization,  we  have 
devised  a  mean  effective  size  (De)  to  represent  the  ice 
crystal  size  distribution  in  radiative  transfer  calculations. 
Based  on  physical  principles,  the  extinction  coefficient 
normalized  by  IWC,  the  single-scattering  albedo,  and  the 
expansion  coefficients  in  the  phase  function  may  be 
expressed  by  a  simple  polynomial  form  in  terms  of  1/De. 
The  coefficients  in  the  polynomials  are  determined  by 
numerical  fitting  to  the  exact  results  obtained  from  the 
light  scattering  and  absorption  programs  using  11 
observed  ice  crystal  size  distributions  (Takano  and  Liou 
1989;  Fu  and  Liou  1993).  The  single-scattering 
parameterization  is  performed  for  six  solar  and  12  IR 
bands.  The  calculations  have  been  carried  out  using  the 
unified  theory  of  light  scattering  by  ice  crystals  developed 
at  the  University  of  Utah  for  all  applicable  size  parameters 
based  on  a  combination  of  bullet  rosettes,  hollow 
columns,  and  plates. 

2.3.  Ice  Microphysks  Parameterizations 

Ice  water  relative  fall  speed  v  is  parameterized  in 
terms  of  pressure  as  follows: 

v(p)=v(p*)(p*/p)^,  (2.5) 

where  p*=400  rttb.  The  reference  value  v(p*)  is  taken  from 
the  parametric  relationship  developed  by  Starr  and  Cox 
(1985). 

Sublimation  of  vapor  to  ice  occurs  when 
qo+q>coq„  where  Co=1.2  when  i=0  andCo=1.05  when 
f  >0.  The  amount  of  water  vapor  sublimated  to  ice  (Aq)  is 
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solved  by  iteration  from  the  following  simultaneous 
equations: 

qo+q-Aq=Coq,(To+r).  (2.6) 

LAq  =  Cp(r-T),  (2.7) 

where  T”  is  the  new  temperature  after  phase  change;  Co  is 
taken  the  same  value  as  given  before.  Then  the 
sublimation  between  water  vapor  and  ice  water  (C)  and  the 
potential  temperature  tendency  due  to  phase  changes  of 
water  (Qc)  can  be  obtained  from 

C  =  Aq/At.  (2.8) 

where  nostpo/p,)"^*^';  pr=1000  mb;  and  At  is  the  time 
interval. 

Sublimation  of  ice  to  vapor  is  diagnosed  when  r>0 
and  qo+q,  <C|q„  where  c,  =  i.05.  The  amount  of  ice 
sublimated  to  water  vapor  is  diagnosed  as  A^ = ,  where  c 
is  a  coefficient  related  to  the  saturation  ratio  and  0  <  c  <  1. 
The  terms  C  and  Qc  can  be  calculated  from  Eqs.  (2.8)  and 
(2.9),  but  Aq  must  be  set  to  be  -Ar.  The  phase  change 
diagnoses  ate  performed  at  five  minutes  intervals  during  a 
simulation. 

2.4.  Boundary  Conditions 

Cyclic  boundary  conditions  are  specified  at  the 
lateral  boundaries.  The  upper  and  lower  boundaries  ate  set 
so  that  the  evolving  disturbances  remain  isolated  within 
the  domain,  i.e., 

;=o.  v=o.  r=o, 
and 

^=_w  i32. 

3t  °  02  ’  3t  °  dz 

For  the  eddy  diilusion  terms,  the  cross-boundary  transport 
is  simply  specified  as  zero  as  follows: 

dz  dz  dz  dz 

The  second-order  moments  are  set  to  zero  at  the  vertical 
boundaries. 

Comprehensive  numerical  schemes  have  been 
designed  for  the  simulation  of  cloud  formation  but  for 
brevity,  they  are  not  given  here. 

3.  CIRRUS  CLOUD  SIMULATION 

Figure  la  shows  the  vertical  profile  of  the 
horizontally  averaged  IWC  at  four  time  steps.  At  t-10 
min,  IWC  increases  rapidly  to  its  maximum  value  due  to 
the  initial  dynamic  and  thermodynamic  forcing.  As  ice 
water  falls  into  the  drier  lower  region,  sublimation  occurs 
and  the  simulated  IWC  decreases.  Toward  the  end  of  the 
simulation,  IWC  reaches  to  a  steady  state.  Also  shown  is 
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the  maximum  value  of  the  simulated  horizontally  averaged 
IWC  which  is  about  0.044  g/m^  at  t=10  min,  much  larger 
than  the  value  0.027  g/m*  reported  by  Starr  and  Cox 
(198S).  The  potential  temperature  tendenccy  due  to  phase 
changes  (Qc)  maintains  positive  in  the  ice  water  region, 
implying  a  net  ice  water  formation  in  this  region.  In  the 
lower  regioa  however,  a  net  ice  destraction  is  found  (Qc  < 
0).  At  first,  formation  is  much  stronger  than  destruction, 
but  at  t=20  min,  destruction  surpasses  generation.  Finally, 
the  generation  of  ice  water  is  balanced  by  the  destraction 
and  the  transport  to  the  lower  level.  The  maximum 
radiaiive  cooling  is  located  at  the  cloud  top,  while  the 
subcloud  region  shows  radiative  wanning.  Radiative 
cooling  and  wanning  decrease  with  time  and  reach  to  a 
quasi-stcady  state  at  the  end  of  the  simulation  (Fig.  Ic). 


Fig.  I.  Vertical  profiles  of  horizontally  averaged  IWC  (a), 
potential  temperature  tendency  due  to  phase  changes  (b), 
and  radiation  (c)  at  four  time  steps. 

During  the  cirrus  cloud  formation  period,  the  latent 
heat  release  in  the  ice  water  region  is  much  larger  than  the 
net  radiative  cooling  in  response  to  the  initial  conditions. 
After  about  t-20  min  of  the  simulation,  these  two  value 
become  comparable,  indicating  that  the  radiative  processes 
begin  to  play  a  significant  role  when  a  sufficient  amount  of 
ice  water  is  produced  (Fig.  lb  and  Fig  Ic). 


Rg.2.  Image  of  IWC  (in  mg/m^)  at  ^10  min. 


Contour  map  of  IWC  at  t=10  min  is  shown  in  Fig. 
2.  Updraft  as  large  as  40  cm/s  is  simulated  and  the 
maximum  simulated  IWC  located  at  about  7  km  is  about 
0.055  g/m\  both  match  the  observed  values  (Heymsfield 
1975) . 


4.  INTERACTIONS  OF  RADUTION, 
TURBULENCE,  AND  ICE  MICROPHYSICS 


IWC(9An**3)  T#(np*i«tut*T*ne«ney(K«iy) 

Fig3.  Vertical  profiles  of  the  horizontally  averaged  IWC 
aat  three  time  steps  (a);  potential  temperature  tendency  due 
to  phase  changes  (Q*)  and  radiation  (Q,)  at  t=40  min 
during  simulations  with  a(solid)  and  without  (dashed) 
turbulence. 

Vertical  profiles  of  the  horizontally  averaged  IWC  at 
three  time  steps  as  well  as  the  latent  heat  and  net  radiative 
heating/cooling  at  t=M0  min  with  and  without  turbulence 
are  presem  in  Fig.  3.  The  primary  role  of  turbulence  is  to 
modulate  the  supersaturation  condition  as  to  affect  ice 
nucleation  ami  consequently  the  rate  at  which  water  vapor 
is  converted  to  ice.  When  the  cirrus  cloud  begins  to  form, 
sublimation  of  water  vapor  to  ice  water  is  primarily 
determined  by  the  initial  forcings.  Turbulent  and  radiative 
processes  appear  to  be  not  important  and  so  the  differences 
between  the  two  experiments  are  small.  After  t=20  min, 
the  direct  effects  of  turbulence  on  phase  changes  become 
apparent.  Qc  increases  by  about  3*C/day  due  to  the 
presence  of  turbulence  and  the  corresponding  horizontally 
averaged  IWC  increases  by  2  mg/m*  (Fig  3a).  At  this 
stage,  however,  the  influence  of  turbulence  on  radiative 
processes  is  still  absent  (figures  are  not  shown  here). 
During  the  maintenance  period  of  the  cirrus  life  cycle, 
feedback  of  turbulence  to  radiative  heating/cooling  becomes 
cignificant  (Fig.  3b).  Due  to  the  including  of  turbulence  in 
the  model,  the  maximum  net  radiative  cooling  in  the  ice 
water  region  is  enhanced  by  about  1  ®C/day,  which  is  25% 
of  the  total  cooling  rate.  At  the  same  time,  Qc  increases  by 
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about  6®C/day.  As  a  result,  IWC  increases  by  about  5 
mg/m^  (Fig.  3a).  Fig.  3b  reveals  that  the  rate  of  net  ice 
formation  at  this  level  is  greater  with  the  tutbulence.  which 
is  consistent  with  the  IWC  difference.  Also,  turbulence 
modifies  the  static  stability  through  radiative  processes,  it 
is  seen  that  the  radiative  stabilization  in  the  upper  portion 
of  the  cloud  generation  region  is  more  apparent  in  the  case 
without  tutbulence. 

We  have  also  carried  out  two  other  experiments 
with  the  inclusion  of  turbulence  but  with  and  without 
incorporating  radiation  in  the  model.  The  results  are 
shown  in  Fig.  4.  Radiation  effects  appear  to  be  not 
substantial  at  first.  However,  at  about  t=20  min  in  the  ice 
water  region,  radiation  enhances  the  latent  heat  by  about 
8*C/day,  implying  that  the  condensation  and  hence  the 
increase  of  the  rate  of  net  formation  are  both  stren^ened 
through  both  direct  effects  from  ice  microphysics  and 
indirect  effects  from  the  induced  entrainment  and  mixing 
(Fig.  4b).  The  corresponding  IWC  differences,  however, 
are  much  smaller  than  that  produced  from  the  turbulent 
effect  (Fig.  4a).  Comparing  Fig.  3a  and  Fig.  4a,  we  see 
that  the  induced  tutbulence,  which  feeds  back  to  the 
radiative  processes,  is  more  significant  in  the  formation  of 
ice  water,  as  evidenced  from  the  larger  IWC  differences. 


Fig.4.  Vertical  proflIes  of  the  horizontally  averaged  IWC 
(a)  and  latent  heat  (b)  for  experiments  with  (solid)  and 
without  (dashed)  radiatioa 


In  suirunaiy,  the  simulation  results  in  terms  of  a 
number  of  physied  parameters  denoted  above  are  sensitive 
to  turbulent  and  the  radiative  processes.  Turbulence 
directly  affects  ice  microphysics  by  producing  higher 
supersaturation.  Also,  it  can  influence  radiation  by 
enhancing  the  net  cooling  rate  and  reducing  the  radiative 
stabilization. 


5.  CONCLUSIONS 

A  two-dimensional  numerical  model  has  been 
developed  to  investigate  the  interactions  among  radiative 


processes,  turbulence,  and  ice  microphysics  and  their 
effects  on  the  formation  of  cirrus  clouds.  Numerical  results 
show  that  cirrus  clouds  are  first  formed  by  the  dynamic  and 
thermodynamic  forcings.  Once  a  sufficient  amount  of  ice 
water  is  produced  in  the  atmosphere,  radiative  processes 
begin  to  play  an  important  role  through  its  net  cooling 
effect  on  the  cloud  top.  The  induced  tutbulence  has 
significant  effects  on  the  phase  changes  and  the  ice  water 
formation.  Also,  tutbulence  affects  the  radiative  processes 
by  increasing  the  net  radiative  cooling  in  the  ice  water 
generation  region,  which  is  important  in  the  maintenance 
of  cirrus  clouds. 
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ABSTRACT 

A  two*dimensional,  time-dependent  numerical  model 
has  been  developed  for  the  studies  of  the  interactions 
among  turbulence,  radiation,  and  microphysics  in  cirrus 
clouds.  The  model  focuses  on  cloud-scale  processes  with  a 
100  m  grid  resolution  and  incorporates  parameterizations 
that  account  for  phase  changes  of  water  and  the  effects  of 
microphysical  structure  on  the  ice  water  vertical  flux  .  Also 
included  is  an  advanced  interactive  radiative  transfer 
scheme,  which  employs  the  6-four-slream  approximation 
for  radiative  transfer,  the  correlated  k-distnbution  method 
for  nongray  gaseous  absorption,  and  the  scattering  and 
absorption  properties  of  hexagonal  ice  crystals.  To 
investigate  the  interactions  of  turbulence,  radiative 
transfer,  and  ice  microphysics  and  their  contributions  to 
the  evolution  of  cirrus  clouds,  a  complete  version  of  the 
second-order  turbulence  closure  model  that  includes 
prognostic  equations  for  all  the  second  moments  has  been 
developed.  Numerical  simulation  of  a  thin  cirrostratus 
shows  that  the  radiative  processes  begin  to  play  an 
important  role  once  a  sufTicient  amount  of  ice  water  is 
produced.  The  induced  turbulence  has  significant  effects  on 
the  phase  changes  and  the  formation  of  ice  water. 
Turbulence  also  enhances  the  radiative  net  cooling  in  the 
ice  water  region,  which  is  an  important  mechanism  for  the 
maintenance  of  cirrus  clouds. 


cirrus  clouds  still  remain  poorly  understood.  Although  it 
is  known  that  radiative  cooling  (and  indirectly  mixing  and 
cloud-top  entrainment)  is  dependent  on  cloud 
microphysics,  the  effect  of  the  induced  turbulence  on  ice 
microphysics,  which  can  feed  back  to  the  radiative  cooling, 
has  not  been  previously  examined.  This  feedback  process 
could  have  an  important  role  in  the  maintenance  of  cirrus 
clouds  (Curry  1985). 

In  this  paper,  we  develop  a  two-dimensional,  time- 
dependent  numerical  model  to  investigate  the  interactions 
among  radiative  processes,  turbulence,  and  ice 
microphysics,  and  their  effects  on  the  formation  and 
maintenance  of  a  cirrus  cloud.  An  advanced  interactive 
radiative  transfer  scheme  for  cirrus  clouds  (Fu  and  Liou 
1993)  has  been  included  in  the  model.  In  addition,  a 
complete  version  of  the  second-order  turbulence  closure 
model  (Level-4  model)  which  contains  the  prognostic 
equations  for  all  of  the  second  moments  has  been 
developed  and  incorporated  into  the  present  model  for 
interaction  and  feedback  studies.  In  Section  2,  we  describe 
the  numerical  model  that  includes  turbulent  transfer  and 
radiative  processes.  Simulation  results  are  presented  in 
Section  3.  In  Section  4,  we  discuss  the  interactions  of 
various  physical  processes  and  their  effects  on  the  life 
cycle  of  a  cirrus  cloud.  Summary  and  conclusions  are  given 
in  Section  5. 

2.  MODEL  DESCRIPTION 


1.  INTRODUCTION 

Ubiquitous  cirrus  clouds  play  an  important  role  in  the 
radiation  field  of  the  carth-almosphcrc  system,  and  hence 
the  earth's  climate.  However,  fundamental  understanding  of 
the  mechanisms  for  their  formation  and  dissipation  is  still 
extremely  limited  and  has  received  substantial  scientific 
attention  of  late  (Liou  1992). 

Several  attempts  have  been  made  to  develop  numerical 
models  to  investigate  the  role  of  various  physical 
processes  in  the  life  cycle  of  a  cirrus  cloud.  The  results  of 
Starr  and  Cox  (1985)  show  that  the  effects  of  radiative 
processes  and  vertical  transports  are  both  significant  in  the 
cirrus  cloud  formation  and  maintenance.  Other  cirrus  cloud 
modeling  studies  (Heymsfield  et  al.  1989;  Jensen  et  al. 
1994)  used  simplified  parameterizations  for  radiative 
processes  and  have  not  considered  turbulence  effects  in  the 
model.  Al  this  point,  interactions  among  radiative 
processes,  turbulence  and  ice  microphysics  and  how  they 
conuibutc  to  the  formation,  maintenance  and  dissipation  of 
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2.1.  Governing  Equations 


The  mean  field  equations  governing  the  motions  in  the 
x-z  plane  based  on  the  Boussinesq  approximation  can  be 
written  as  follows: 

at  m  Uq 
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where  the  subscripts  refer  to  the  initial  state,  which  is 
horizontally  unifonn,  cloud  free  and  at  rest  except  for  a 
constant  vertical  motion  Wq;  the  unsubscripted  variables 
are  regarded  as  mean  Helds  based  on  the  initial  state; 
variables  with  prime  represent  turbulent  fluctuations;  Qr 
and  Qc  are  the  potential  temperature  tendencies  due  to 
radiative  processes  and  phase  changes  of  water, 
respectively;  v  is  the  ice  water  relative  fall  speed;  C  is  the 
sublimation  between  water  vapor  and  ice  water,  q  is  water 
vapor  specific  humidity;  t  is  ice  water  mixing  ratio;  v  is 
eddy  viscosity;  ke,  k,,  k,  arc  thermal,  water  vapor  and 

ice  water  eddy  diflusivities,  respectively;  £’  is  related  to 
the  ratio  of  molecular  weight  of  water  to  that  of  dry  air,  and 
other  notations  are  conventional. 


2.2.  Turbulent  Flux  Equations 


The  governing  equations  (2.1)-(2.4)  require  nine 
turbulent  fluxes,  which  can  be  predicted  using  the  second- 
order  closure  based  on  12  prognostic  equations  for  the 
second  moments.  The  closure  assumptions  are  described  in 
Mellor  and  Yamada  (1974).  The  mean  Reynolds  stress 
model  equations  are  given  by 
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Equation  (2.5)  contains  three  prognostic  equations  for 
u*^,  u’w*,  and  w*^,  while  Eq.  (2.7)  contains  two 
equations  for  u’0'  and  w‘0\  The  turbulent  fluxes  with 
respect  to  q  and  f  can  be  expressed  by  following  Eqs.  (2.6) 
and  (2.7)  and  arc  not  given  here.  In  these  equations,  P  is 
the  coeflicient  of  thennal  expansion;  gjK0,0,-g)  the  gravity 
vector,  q*  is  the  square  root  of  the  turbulence  kinetic 
energy;  c/  is  a  constant  taken  from  Mellor  and  Yamada 
(1974);  and  X„  Xj,  Xj,  f,  fj,  and  a,  Aj  are  all  length 
parameters  which  are  related  to  length  scale  1. 

For  the  length  scale,  we  have  adopted  the 
interpolation  formula  developed  by  Blackdar  (1962)  in  the 
fomfi 


1  +  kz/lo 

where  k  is  von  Karman  constant;  Iq  =  0.00027  *  Uq  /f ;  f  is 


the  Coriolis  parameter,  and  Uq  is  the  geostrophic  wind 
velocity,  which  is  set  at  a  constant  value  16  m/s. 


We  have  used  the  radiation  parameterization  scheme 
developed  at  the  University  of  Utah  for  the  cirrus  cloud 
modeling  studies.  This  scheme  is  based  on  the  delta-four- 
stream  approximation  developed  by  Liou  ct  al.  (1988).  It 
can  be  applied  to  a  nonhomogeneous  cloudy  atmosphere 
by  dividing  it  into  a  number  of  layers  each  of  which  can  be 
considered  homogeneous  with  respect  to  the  single¬ 
scattering  albedo  and  phase  function.  Using  the  boundary 
conditions  at  the  top  and  bottom  of  the  atmosphere  and  the 
continuity  requirements  at  the  interface  of  the  predividrf 
homogeneous  layers,  the  unknown  coefficients  in  the 
analytic  solution  for  the  transfer  equation  can  be 
determined,  leading  to  the  calculations  of  fluxes. 

A  parameterization  of  non-gray  gaseous  absorption 
based  on  the  correlated  k-distribution  method  has  been 
developed  by  Fu  and  Liou  (1992).  In  this  method,  the 
cumulative  probability  g  of  the  absorption  coefficient  in  a 
^ctral  interval  is  used  to  replace  wavenumber  as  an 
independent  variable.  Monochromatic  solutions  of  the 
radiative  transfer  equation  for  the  multiple-scattering 
atmosphere  can  be  used  in  connection  with  a  given  g.  The 
minimum  number  of  quadrature  points  in  the  g  ^ce  ranges 
from  one  to  about  ten  for  different  absorbing  gases  in 
different  spectral  regions. 

The  solar  and  IR  spectra  are  divided  into  a  number  of 
spectral  intervals  according  to  the  location  of  the  gaseous 
absorption  bands.  Absorption  due  to  H2O,  O3,  CO2,  and 
O2  is  accounted  for  in  the  solar  spectrum.  In  the  thermal 
infrared,  absorption  due  to  H2O,  O3,  CO2,  CH4,  and  N2O  is 
included.  Using  the  correlated  k-distribution  method,  121 
spectral  calculations  arc  required  for  each  vertical  profile. 

A  new  parameterization  has  been  developed  for  the 
calculations  of  the  single-scattering  properties  for 
hexagonal  ice  crystals.  In  this  parameterization,  we  have 
devised  a  mean  effective  size  (De),  which  is  defined  as  the 
mean  width  weighted  by  the  cross-section  area,  to 
represent  the  ice  crystal  size  distribution  in  radiative 
transfer  calculations.  Based  on  physical  principles,  the 
extinction  coefficient  normalized  by  IWC,  the  single- 
scattering  albedo,  and  the  expansion  coefficients  in  the 
phase  function  may  be  expressed  by  a  simple  polynomial 
form  in  terms  of  l/D^.  The  coefficients  in  the  polynomials 
are  determined  by  numerical  fitting  to  the  exact  results 
obtained  from  the  light  scattering  and  absorption  programs 
using  11  observed  ice  crystal  size  distributions  (Takano 
and  Liou  1989;  Fu  and  Liou  1993).  The  single-scattering 
parameterization  is  performed  for  six  solar  and  12  IR  bands. 
The  calculations  have  been  carried  out  using  the  unified 
theory  of  light  scattering  by  ice  crystals  developed  at  the 
University  of  Utah  for  all  applicable  size  parameters  based 
on  a  combination  of  bullet  rosettes,  hollow  columns,  and 
plates. 

The  radiation  parametcrizaion  program  is  driven  by 
the  mean  effective  size  and  ice  water  path,  which  is  the 
product  of  IWC  and  the  cloud  thickness.  In  the  present 
cirrus  cloud  model,  the  ice  crystal  size  is  not  predicted  and 
is  prescribed  using  a  typical  cirrostratus  value. 

2.4.  Ice  Microphysics  Parameterizations 

Ice  water  relative  fall  speed  v  is  parameterized  in  terms 
of  pressure  as  follows: 

v(p)  =  v(p*)  (p*/p)^. 


2.3.  Radiative  Transfer  Parameterizations 


(2.9) 


where  p^=400  mb.  The  reference  value  v(p^)  is  taken  from 
the  parametric  relationship  developed  by  Starr  and  Cox 
(1985). 

Sublimation  of  vapor  to  ice  occurs  when  qo  +<)  >  Coqg* 

where  Co=1.2  when  f=0  and  Co==1.05  when  i>0.  The 
amount  of  water  vapor  sublimated  to  ice  (Aq)  is  solved  by 
iteration  from  the  following  simultaneous  equations: 

qo+q-Aq=Coq,(To-i-T*).  (2.10) 

LAq=c^(r-T).  (2.11) 

where  T”  is  the  new  temperature  after  phase  change;  c©  is 
taken  the  same  value  as  given  before.  Then  the  sublimation 
between  water  vapor  and  ice  water  (C)  and  the  potential 
temperature  tendency  due  to  phase  changes  of  water  (Qc) 
can  be  obtained  Grom 

C  =  Aq/Al,  (2.12) 

(2.3) 

where  Hq  =(Po/Pf)^^^;  Pr=1000  mb;  and  At  is  the  time 
interval. 

Sublimation  of  ice  to  vapor  is  diagnosed  when  />0 
and  qo’i’qi  <<^iqit  where  C|  =:1.05.  The  amount  of  ice 
sublimated  to  water  vapor  is  diagnosed  as  A^  =  cf ,  where  c 
is  a  coefticient  related  to  the  saturation  ratio  and  0  ^  c  ^  1. 
The  terms  C  and  Qc  can  be  calculated  from  Eqs.  (2.12)  and 
(2.13),  but  Aq  must  be  set  to  be  -Af .  The  phase  change 
diagnoses  arc  performed  at  five  minutes  intervals  during  a 
simulation. 


2.S.  Boundary  Conditions  And  Numerical  Methods 


Cyclic  boundary  conditions  are  specified  at  the  lateral 
bound^es.  The  upper  and  lower  boundaries  are  set  so  that 
the  evolving  disturbances  remain  isolated  within  the 
domain,  i.e., 

5=0,  ysO,  f=0. 
and 


ii--w  ii-w  isi 

dt  °  dz  '  9t  ®  3z 

For  the  eddy  diffusion  terms,  the  cross-boundary  transport 
is  simply  specifled  as  zero  as  follows: 


dz  dz  dz  dz 

The  second-order  moments  are  set  to  zero  at  the  vertical 
boundaries. 

The  dimensions  of  the  computational  domain  used  are 
3.2  km  in  the  vertical  and  6.4  to  in  the  horizontal  with  a 
uniform  horizontal  and  vertical  resolution  of  100  m.  The 
nonlinear  term  of  the  vorticity  equation  (2.1)  is  evaluated 
by  the  finite-difference  scheme  in  terms  of  the  (lux  form 
using  Arakawa’s  (1966)  nine-point  Jacobian,  which 
conserves  both  the  mean  kinetic  energy  and  mean  square 
vorticity.  The  nonlinear  terms  in  equations  (2.2)  and  (2.3) 
are  computed  following  the  method  developed  by  Lilly 
(1965).  The  terms  in  Eq.  (2.4)  are  obtained  by  an  upstream 
flux  scheme.  Equations  (2.1H2.3)  are  solved  by  using  the 
second-order  Adams-Bashfoi^  scheme,  except  that  the 
Euler  forward  scheme  is  used  for  the  initial  time  step.  The 
Euler  forward  scheme  is  also  used  for  the  solution  of  ice 
water  denoted  in  Eq.  (2.4)  to  ensure  computational 
stability. 


3.  CIRRUS  CLOUD  SIMULATION 

The  initial  thermal  stratifications  used  in  the  present 
model  are  a  number  of  specified  lapse  rates  in  various 
vertical  regions,  which  are  expressed  in  terms  of  the 
pseudo-adiabatic  lapse  rate  governing  the  unmixed  vertical 
displacements  at  ice  saturation,  typical  of  thin  cirrus  cloud 
layers.  There  is  no  ice  water  initially.  The  disturbance  is 
initiated  through  the  prescribed  potential  temperature 
perturbation  at  all  grid  points  in  the  layer  from  z=^.55  to 
7.55  to,  with  a  maximum  and  a  mean  value  of  O.l  and 
0.05*C,  respectively. 


Fig.l.  Vertical  profiles  of  horizontally  averaged  IWC  (a), 
the  potential  temperature  tendency  due  to  phase  changes 
(b),  and  radiation  (c)  at  four  time  steps. 

Figure  la  shows  the  vertical  profile  of  the 
horizontally  averaged  IWC  at  four  time  steps  during  a 
cirrus  cloud  simulation.  At  t^lO  min,  IWC  increases 
rapidly  to  its  maximum  value  due  to  the  initial  dynamic  and 
thermodynamic  forcings.  Because  the  values  of  the  ice 
water  relative  fall  spe^  at  that  time  are  larger  than  the 
updraft  speeds,  ice  water  falls  to  levels  below  and  the 
location  of  the  maximum  IWC  shifts  from  7.15  to  at  t=10 
min  to  6.95  to  at  t=20  min.  As  ice  water  falls  into  the  drier 
lower  region,  sublimation  occurs  and  the  simulated  IWC 
decreases.  Toward  the  end  of  the  simulation,  IWC  reaches 
to  a  steady  state  in  which  the  formation  and  destruction  of 
ice  water  is  balanced.  Also  shown  is  the  maximum  value  of 
the  simulated  horizontally  averaged  IWC  which  is  about 
0.044  g/m^  at  t=10  min,  located  at  about  7  km,  much  larger 
than  the  value  0.027  g/m^  reported  by  Starr  and  Cox 
(1985).  These  processes  are  also  clearly  illustrated  from  the 
vertical  profile  of  potential  temperature  tendency  due  to 
phase  changes  of  water  (Qc)  (Fig.  lb).  Qc  maintains 
positive  in  the  ice  water  region,  implying  a  net  ice  water 
formation  in  this  region.  In  the  lower  region,  however,  a 
net  ice  destruction  is  found  (Qc  <  0).  At  first,  formation  is 
much  stronger  than  destruction,  but  at  1=20  min, 
destruction  surpasses  generation.  Finally,  the  generation 
of  ice  water  is  balanced  by  the  destruction  and  the 
transport  to  the  lower  level.  The  vertical  profiles  of  net 
radiative  coolingAvarming  at  different  time  steps  are  shown 
in  Fig.  Ic.  The  maximum  radiative  cooling  is  located  at  the 
cloud  top,  while  the  subcloud  region  shows  radiative 
warming.  Radiative  cooling  and  warming  decrease  with 
time  and  reach  to  a  quasi-steady  state  at  the  end  of  the 
simulation. 

Vertical  profiles  of  the  potential  temperature  change 
due  to  phase  change,  solar,  infrared  and  net  radiation  are 
illustrate  in  Fig.  2.  During  the  cirrus  cloud  formation 


period,  the  latent  heat  release  in  the  ice  water  region  is 
much  larger  than  the  net  radiative  cooling  in  response  to 
the  initial  conditions.  After  about  t=20  min  of  the 
simulation,  these  two  value  become  comparable,  indicating 
that  the  radiative  processes  begin  to  play  a  significant  role 
when  a  sufficient  amount  of  ice  water  is  produced. 
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Fig.2.  Vertical  profiles  of  temperature  change  due  to  (a) 
phase  changes,  (b)  solar,  (c)  infrared,  and  (d)  net  radiative 
processes  at  four  time  steps. 

Contour  map  of  IWC  at  t=10  min  is  shoN^n  in  Fig.  3. 
Updraft  as  large  as  40  cm/s  is  simulated  and  the  maximum 
simulated  IWC  located  at  about  7  km  is  about  0.05S  g/m^, 
both  match  the  observed  values  (Heymsfield  1975) . 


Fig.3.  Image  of  IWC  (in  mgAn^)  at  t=10  min. 

From  the  preceding  discussion,  it  is  clear  that  the 
radiative  processes  play  an  important  role  in  the  life  cycle 
of  a  cirrus  cloud,  especially  in  the  maintenance  period.  The 
present  model  also  produces  more  ice  water  than  that 
reported  in  Starr  and  Cox  (1985)  in  whose  model  the  effects 
of  turbulence  and  radiation  have  not  been  treated 
comprehensively. 

4.  INTERACTIONS  OF  RADIATION,  TURBULENCE. 

AND  ICE  MICROPHYSICS 

In  this  section,  we  investigate  the  interactions  and  the 
feedback  of  the  entrainment  and  mixing  to  the  radiative 
cooling  during  the  simulation  of  the  life  cycle  of  a  cirrus 


cloud.  Vertical  profiles  of  the  horizontally  averaged  IWC  at 
three  time  steps  as  well  as  the  latent  heat  and  net  radiative 
heating/cooling  at  t=40  min  with  and  without  turbulence 
are  present  in  Fig.  4.  The  primary  role  of  turbulence  is  to 
modulate  the  supersaturation  condition  as  to  affect  ice 
nucleation  and  consequently  the  rate  at  which  water  vapor 
is  converted  to  ice.  \^en  the  cirrus  cloud  begins  to  form, 
sublimation  of  water  vapor  to  ice  water  is  primarily 
determined  by  the  initial  forcings.  Turbulent  and  radiative 
processes  appear  to  be  not  important  and  so  the  differences 
between  the  two  experiments  are  small.  After  t=20  min,  the 
direct  effects  of  turbulence  on  phase  changes  become 
apparent.  Qc  increases  by  about  3  ®C/day  due  to  the 
presence  of  turbulence  and  the  corresponding  horizontally 
averaged  IWC  increases  by  2  mg/m^  (Fig.  4a).  At  this  stage, 
however,  the  influence  of  turbulence  on  radiative  processes 
is  still  absent  (figures  are  not  shown  here).  During  the 
maintenance  period  of  the  cirrus  life  cycle  at  t=40  min, 
feedback  of  turbulence  to  radiative  heating/cooling 
becomes  significant  (Fig.  4b).  Solar  heating  is  less  than  IR 
cooling  throughout  all  layers,  resulting  in  net  cooling. 
Due  to  the  including  of  turbulence  in  the  model,  the 
maximum  net  radiative  cooling  in  the  ice  water  region  is 
enhanced  by  about  1  •C/day,  which  is  25%  of  the  total 
cooling  rate.  At  the  same  time,  Qc  increases  by  about 
6®C/day.  As  a  result,  IWC  increases  by  about  5  mg/m^ 
(Fig.  4a).  Fig.  4b  reveals  that  the  rate  of  net  ice  formation  at 
this  level  is  greater  with  the  turbulence,  which  is 
consistent  with  the  IWC  difference.  Also,  turbulence 
modifies  the  static  stability  through  radiative  processes.  It 
is  seen  that  the  radiative  stabilization  in  the  upper  portion 
of  the  cloud  generation  region  is  more  apparent  in  the  case 
without  turbulence. 


IWC  Tamperatufa  Taodeocy  (K/day) 

Fig.4.  Vertical  profiles  of  the  horizontally  averaged  IWC  at 
three  time  steps  (a);  potential  temperature  tendency  due  to 
phase  changes  (QJ  and  radiation  (Q^)  at  t=40  min  during 
simulations  with  (solid)  and  without  (dashed)  turbulence. 

We  have  also  carried  out  two  other  experiments  with 
the  inclusion  of  turbulence  but  with  and  without 
incorporating  radiation  in  the  model.  The  results  are 
shown  in  Fig.  5.  Radiation  effects  appear  to  be  not 
substantial  at  first.  However,  at  about  t=20  min  in  the  ice 
water  region,  radiation  enhances  the  latent  heat  by  about 
8*C/day,  implying  that  the  condensation  and  hence  the 
increase  of  the  rate  of  net  formation  are  both  strengthened 
through  both  direct  effects  from  ice  microphysics  and 
indirect  effects  from  the  induced  entrainment  and  mixing 


(Fig.  5b).  The  corresponding  IWC  differences,  however,  are 
much  smaller  than  that  produced  from  the  turbulent  effect 
(Fig.  5a).  Comparing  Fig.  4a  and  Fig.  5a,  we  see  that  the 
induced  turbulence,  which  feeds  back  to  the  radiative 
processes,  is  more  significant  in  the  formation  of  ice  water, 
as  evidenced  from  the  larger  IWC  differences. 
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Fig.5.  Vertical  profiles  of  the  horizontally  averaged  IWC 
(a)  and  latent  heat  (b)  for  experiments  with  (solid)  and 
without  (dashed)  radiation. 

Another  two  experiments  have  been  made  under 
conditions  without  the  inclusion  of  turbulence,  but  with 
and  without  the  incorporation  of  radiation.  Only  very 
small  differences  in  both  IWC  and  latent  heat  are  found  in 
these  two  experiments.  This  is  due  to  the  fact  that  radiation 
itself  cannot  play  a  significant  role  in  the  net  ice  formation 
in  the  ice  generation  region  nor  in  the  net  sublimation  of 
ice  to  water  vapor  in  the  lower  region.  Likewise,  if 
radiation  is  not  included,  turbulence  cannot  produce  larger 
differences  in  IWC,  as  shown  in  Fig.  4. 

In  summary,  the  simulation  results  in  terms  of  a  number 
of  physical  parameters  denoted  above  are  sensitive  to 
turbulent  and  the  radiative  processes.  Turbulence  directly 
affects  ice  microphysics  by  producing  higher 
supersaturation.  Also,  it  can  influence  radiation  by 
enhancing  the  net  cooling  rate  and  reducing  the  radiative 
stabilization.  The  interactions  among  radiation,  turbulence 
and  microphysics  play  an  important  role  in  the  life  cycle  of 
cirrus  clouds. 

5.  CONCLUSIONS 

A  two-dimensional  numerical  model  has  been 
developed  to  investigate  the  interactions  among  radiative 
processes,  turbulence,  and  ice  microphysics  and  their 
effects  on  the  formation  of  cirrus  clouds.  A  second-order  4- 
level  turbulence  closure  model  has  been  developed  and 
incorporated  into  the  present  model.  An  advanced 
radiative  scheme  is  used  to  account  for  the  radiative 
processes.  Sublimation  between  water  vapor  and  ice  water 
as  well  as  the  fall  spe^  of  ice  water  have  been  accounted 
for  based  on  parametcrizations. 

Gurus  clouds  are  first  formed  by  the  dynamic  and 
thermodynamic  forcings.  Once  a  sufficient  amount  of  ice 
water  is  produced  in  the  atmosphere,  radiative  processes 
begin  to  play  an  important  role  through  its  net  cooling 
effect  on  the  cloud  top.  The  induced  turbulence  has 
significant  effects  on  the  phase  changes  and  the  ice  water 


formation.  Also,  turbulence  affects  the  radiative  processes 
by  increasing  the  net  radiative  cooling  in  the  ice  water 
generation  region,  which  is  important  in  the  maintenance 
of  cirrus  clouds. 
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1.  INTRODUCTION 

In  a  period  of  history  when  jet  aircraft 
dominate  the  skies,  questions  arise  as  to  the  impact  of 
such  activity  on  the  atmosphere.  Machta  and 
Carpenter  (1971)  and  Liou  et  a/  (19901.  showed  that 
there  was  an  increase  in  the  amount  of  high  cloud 
cover  between  1948  to  1984  when  lower  clouds  were 
not  present  at  various  sites  over  the  United  States. 
They  further  suggested  that  a  correlation  existed 
between  this  increase  in  high  cloud  cover  and 
increases  in  jet  traffic.  More  recently,  the  formation 
and  effect  of  contrail  clouds  arising  from  jet  engine 
exhaust  plumes  has  received  considerable  attention. 
Contrail  clouds  trap  outgoing  longwave  radiation 
emitted  from  the  earth  and  the  troposphere  and  keep 
the  surface  warm,  referred  to  as  the  greenhouse  effect. 
At  the  same  time,  these  clouds  reflect  incoming  solar 
radiation  and  reduce  the  solar  energy  available  to  the 
surface,  referred  to  as  the  solar  albedo  effect.  The 
balance  between  these  two  effects  determines  the 
tendency  of  the  surface  temperature  change.  Various 
studies  have  shown  that  cirrus  and  contrail  cirrus  tend 
to  have  ir  greenhouse  effect  dominating  the  solar 
albedo  effect,  therefore  causing  a  net  warming  at  the 
surface  (Manabe  1975,  Freeman  and  Liou  1979,  Sinha 
and  Shine  1 974,  Liou  1 986). 

A  significant  increase  in  jet  fuel  consumption 
in  the  1 960's  corresponding  to  an  increase  in  cirrus 
cloud  cover  has  been  documented.  While  for  the 
current  study  an  accurate  figure  of  jet  fuel 
consumption  over  the  local  areas  of  cloud  data  analysis 
is  not  available,  extensive  studies  on  global  jet  fuel 
consumption  and  air  traffic  patterns  and  populations 
have  been  undertaken.  Jet  fuel  consumption  statistics 
and  projections  between  1978-2000  showing  continual 
increases  Nusser  and  Schmitt  1990),  which  implies 
increased  tropospheric  water  vapor  emissions.  While 
contrail  formation  is  a  subject  requiring  considerable 
research,  it  is  clearly  related  to  water  vapor 
concentration.  This  could  be  one  reason  that  as  jet 
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traffic  appears  to  be  increasing,  so  does  cirrus  cloud 
extent.  Consequently,  a  correlation  between  jet  air 
traffic  arvl  cirrus  cloud  cover  becomes  a  tool  in 
assessing  the  climatic  impact  of  the  former. 

Contrails  tend  to  form  in  the  vicinity  of  the 
tropopause.  In  this  is  the  region  the  ambient 
temperature  allows  supersaturated  water  vapor 
emitted  from  jet  engines  to  quickly  go  into  ice  phase. 
Residency  times  for  particles  may  vary  from  days 
below  the  tropopause  up  to  a  year  above  the 
tropopause  in  the  lower  stratosphere  (GrassI  1990). 
As  the  ambient  water  vapor  mixing  ratio  is  low  in  the 
upper  troposphere  and  lower  stratosphere,  a  small 
amount  of  water  vapor  in  the  exhaust  plume  has  a 
substantially  greater  impact  at  these  altitudes  than  in 
the  lower  troposphere.  Thus  jet  aircraft  could  be  a 
major  source  of  water  vapor  in  this  region  of  the 
atmosphere.  Further  analysis  of  air  traffic  patterns 
shows  that  in  general  cruising  altitudes  are  within 
contrail  cirrus  producing  levels.  There  are  other 
considerations  as  well,  the  emission  of  SO.  and  NO, 
from  jet  engines  could  contribute  to  ice  nucleation  as 
well  as  have  their  own  radiative  effects.  And  even  if 
the  water  vapor  emitted  does  not  form  contrail  cirrus 
it  is  the  most'  important  feedback  process  involved  in 
greenhouse  perturbations.  Water  vapor  emission  from 
jet  engines  is  roughly  proportional  >to  jet  fuel 
consumption!  Nusser  and  Schmitt  1990),  and 
projections  indicate  that  global  jet  fuel  consumption  will 
continue  to  increase  into  the  twenty  first  century. 

2.  CIRRUS  CLOUD  CLIMATOLOGY 

While  cirrus  cloud  data  is  not  prevalent,  this 
paper  studies  some  of  the  existing  climatologies. 
Cirrus  cloud  cover  data  was  analyzed  for  the  later 
twentieth  century.  The  data  is  comprised  of  local 
surface  observations  over  urban  airports  in  the  United 
States  spanning  forty  years,  and  five  years  of  satellite 
data  taken  from  the  HIRS  radiometer  on  the  NOAA 
polar  orbiting  satellite. 

The  surface  observations  of  cloud  covers 
include  datasets  for  Salt  Lake  City,  Chicago,  Denver, 
and  San  Francisco.  To  be  consistent  with  previous 
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cloud  cover  data  sets  (Machta  and  Carpenter  1971, 
Liou  et  al.  1990),  cirrus  data  was  only  considered 
when  the  combined  cloud  cover  of  lower  and  middle 
clouds  was  less  than  or  equal  to  5%.  This  was  held 
for  all  the  surface  observation  data  sets.  Observations 
were  taken  hourly  and  from  this  data  a  daily  average, 
monthly,  and  annual  average  cloud  covers  were 
calculated.  Subject  to  the  preceding  conditions,  all 
observed  layers  containing  the  cirrus  types  described 
below  were  summed  to  obtain  the  total  cirrus  cloud 
cover  for  that  observation.  It  should  be  noted  that  this 
analysis  would  tend  to  underestimate  cirrus  covers  due 
to  overlap  over  underlying  opaque  low  and  middle 
cloud  layers.  The  cirrus  cloud  cover  statistics 
presented  here  considers  various  cloud  types:  cirrus, 
cirrostratus,  cirrocumulus,  and  cirrocumulus  standing 
lenticular. 

Based  on  the  analysis  of  domestic  jet  fuel 
consumption,  significant  increases  occurred  in  about 
1964  (Liou  et  al.  1990).  Thus  we  study  the  cloud 
climatology  for  various  midlatitude  cities  using  this  year 
as  a  reference  year.  Viewing  the  Salt  Lake  City 
observations,  two  populations  of  cloud  data  can  be 
discerned,  as  illustrated  in  Rgure  1(a).  The  first 
covers  the  dates  between  1 949*1 964  and  has  an 
average  values  of  0.115.  Then  the  cloud  cover 
increases  from  1964-1 993  with  an  average  of  0.202. 
A  linear  fit  is  also  made  for  the  whole  record  to  show 
the  tendency.  In  effect  we  see  a  doubling  of  cloud 
cover  over  the  time  period  of  the  entire  data  set.  While 
cloud  cover  is  variant  under  weather  and  temporary 
climate  perturbations  such  as  volcanic  aerosol,  the  data 
shows  a  significant  increase  in  annual  average  cloud 
cover  as  is  indicated  in  the  figure. 

Another  dataset  of  comparable  time  period 
was  investigated.  Chicago  shows  increased  annual 
average  values  of  cirrus  cloud  cover  ( Fig.  1  (b)).  though 
not  as  strong  an  increase  as  Salt  Lake  City  perhaps  due 
to  the  more  saturated  environment  of  the  Chicago 
area..  The  data  set  for  Denver  dates  back  to  1 945 
but  no  data  is  available  for  the  1970*1 984  period. 
Significant  differences  in  mean  cloud  cover  are  evident 
before  and  after  about  1964,  as  shown  in  Fig  1  (c).  For 
San  Francisco,  data  is  available  after  1985  and  the 
cirrus  cloud  cover  values  are  comparable  to  those  for 
Salt  Lake  City  and  Denver,  it  should  be  noted  that 
within  a  fifty  year  period  the  data  is  higNy  variable  for 
all  the  cases  and  hence  it  is  appropriate  to  examine  the 
cloud  cover  tendencies. 

Supplementing  the  surface  observation  cirrus 
cloud  cover  are  datasets  determined  from  the  polar- 
orbiting  High-Resolution  Infrared  Radiation  Sounder 
(HIRS)  on  board  NOAA  satellites,  over  the  time  period 
from  June  1989  to  September  1995.  The  HIRS 
observations  are  based  on  the  so  called  C02  slicing 
method  which  uses  the  multispectral  channels  in  the 
15  mm  CO2  band  to  retrieve  cloud  cover  statistics 
(Wylie  et  al.  1 994).  Due  to  the  method  of  detection, 
the  cirrus  classification  is  inherently  different  between 
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Figure  1 .  Cloud  covers  for  respective  cities. 

surface  and  satellite  observations.  The  HIRS  data  is 
expressed  in  terms  of  effective  emissivity,  which  is  an 
average  emissivity  measured  over  an  approximately  20 
km  X  20  km  field  of  view.  For  the  HIRS  observations 
to  be  considered  cirrus  the  effective  emissivity  must 
be  less  than  0.95,  corresponding  to  an  optical  depth 
less  than  3.0.  HIRS  data  shows  cirrus  clouds 
regardless  of  lower  clouds.  This  will  tend  to  make  the 
HIRS  cloud  cover  statistics  have  higher  values  than  the 
surface  observation  statistics.  The  portion  of  HIRS 
data  included  in  the  present  study  covers  the  65*S  to 
SS^N  geographical  area.  Poleward  of  this  band  was 
not  considered,  due  to  data  availability.  This  dataset 


Figure  2.  Cirrus  cloud  covers  detected  by  HIRS. 

includes  cirrus  in  two  categories;  thin  cirrus  with 
effective  emissivity  of  less  than  or  equal  to  0.5 
(optical  depth  of  0.7  or  smaller),  and  thick  cirrus  with 
effective  emissivity  of  greater  than  0.5  (optical  depth 
greater  than  0.7).  The  data  was  presented  in  the  form 
of  monthly  averages  in  which  measurements  made 
over  water  and  land  were  distinguished.  For  the 
purposes  of  this  study  the  measurements  over  water 
were  weighted  75%  against  the  25%  over  land  to  form 
a  single  average  cirrus  cover  for  each  month. 

Figure  2  shows  the  total,  thick  and  thin  cirrus 
covers  detected  by  HIRS.  In  calculating  the  total  cirrus 
cloud  cover,  the  thin  cirrus  frequency  was  multiplied  by 
0.72,  a  statistical  average  of  the  effective  emissivity, 
to  obtain  the  cloud  fraction.  While  significantly  higher 
than  surface  observation  data  as  discussed  earlier,  the 
HIRS  data  shows  a  global  increasing  trend  over  the 
four  year  period  of  June  1 939  to  September  1993.  In 
June  1989  the  value  is  0.30.  The  cirrus  cover 
increases  over  time  to  finish  at  an  average  of  0.40. 
Hence  local  surface  observations  taken  over  a  long 
time  scale  and  satellite  observations  taken  over  a  large 
spatial  scale  but  short  time  scale  demonstrate  a 
definite  increase  in  cirrus  cloud  cover  based  on 
consistent  satellite  observations. 

3.  CLIMATIC  IMPACT  OF  CIRRUS  ON 
SURFACE  TEMPERATURE 

To  consider  the  effects  of  global  cirrus  cover 
to  climate,  a  one-dimensional  radiative-convective 
model  {Frankel  1 995)  was  used  to  calculate 
temperature  perturbations  under  a  condition  of 
changing  cirrus  cover  over  time.  This  model  is 
designed  to  maintain  global  energy  balance  and 
progress  to  a  state  of  thermal  atmospheric  equilibrium 
in  the  context  of  an  annual  and  spatial  averaged  entity. 
This  affords  the  opportunity  to  look  at  atmospheric 
cloud-radiation  interactions  in  great  detail  and  with  a 


Figure  3.  Surface  temperatures  from  model  output. 

high  vertical  resolution.  The  model  uses  the  radiative 
scheme  of  Fu  and  Liou  (1993).  The  scheme  consists  of 
the  delta-four  stream  method  for  the  transfer  of  both 
solar  and  thermal  infrared  radiation  in  inhomogeneous 
atmospheres  including  clouds  and  the  correlated  k- 
distribution  method  to  handle  non-gray  gaseous 
absorption  in  multiple  scattering  atmospheres.  In 
particular,  the  scattering  and  absorption  properties  of 
hexagonal  ice  crystals  were  incorporated  in  the 
radiation  scheme  in  terms  of  parameterization  forms 
based  on  the  unified  theory  of  light  scattering  by  ice 
crystals  covering  all  sizes  and  shapes  (Liou  and  Takano 
1994,  Yang  and  Liou  1996  a,  b).  The  convective  part 
of  the  model  computes  a  convective  adjustment  and 
implicit  convective  flux  similar  to  that  found  in 
traditional  radiative-convective  models  (Manabe  1975). 
The  model  was  tuned  to  the  Standard  Atmosphere  with 
the  total  cirrus  cloud  cover  fixed  to  a  value 
corresponding  to  the  first  month  of  the  HIRS  data. 

To  asses  the  effects  of  the  increase  in  total 
cirrus  cover  the  1  -D  model  was  run  interactively  using 
the  HIRS  results.  Initialized  to  thermal  equilibrium  and 
a  surface  temperature  of  288.07  K  under  the  condition 
of  cirrus  cloud  cover  fixed  to  the  value  from  June  1 989 
HIRS  data,  the  model  was  then  let  to  run  for  a  number 
of  timesteps  equivalent  to  four  years,  updating  the* 
cloud  cover  in  the  model  every  month  to  a  new  HIRS 
value.  The  model  output  temperatures  are  shown  in 
Fig.  3.  A  linear  regression  of  this  curve  shows  an 
increase  of  1 .2  K  for  a  five  year  period.  Liou,  et  ai 
(1990)  found  for  a  5%  increase  in  contrail  cirrus  cover 
rendered  a  surface  temperature  increase  of  1  K.  The 
net  cirrus  increase  in  the  HIRS  dataset  is  approximately 
10%  fractional  cloud  cover,  increasing  from 
approximately  30%  cover  to  40%  cover.  The  present 
result  of  1.2  K  temperature  change  appears  to  be 
consistent  with  other  climate  models.  This  is  a 
substantial  temperature  perturbation,  and  suggests  that 
without  other  influences,  the  increase  of  cirrus  cloud 
cover  could  lead  to  produce  a  significant  warming 
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effect. 


Although  the  conditions  under  which  jet 
aircraft  contrails  become  cirrus  clouds  are  not  well 
understood  at  this  point,  both  surface  and  satellite 
observations  show  that  cirrus  cloud  cover  appears  to 
be  on  the  increase.  Thin  cirrus  clouds  have  a  warming 
effect  upon  the  lower  troposphere  through  radiative 
forcing,  and  a  numerical  model  study  suggests  that 
observable  cirrus  cloud  cover  increase  would  have  a 
noticeable  impact  upon  surface  temperatures. 
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On  the  radiative  properties  of  contrail  cirrus 


K.  N.  Liou‘,  P.  Yang*,  Y.  Takano*,  K.  Sassen*,  T.  Charlock®,  and  W.  Arnott^ 


Abstract*  Using  the  observed  ice  crystal  size  distribu¬ 
tion  in  contrail  cirrua  from  SUCCESS,  we  have  catri^ 
out  the  scattering  and  absorption  calculations  based  on 
a  unified  theory  for  light  scattering  by  ice  crystals  cov. 
ering  all  siza  and  shapes.  We  illustrate  the  effects  of 
ice  crystal  size  and  surface  roughness  on  the  ecatter- 
phase  function  features  for  remote  sensing  applicsr 
tions-  The  extinction  coefficient  and  single-scattering 
albedo  exhibit  a  minimum  feature  at  2.85  referred 
to  as  the  Christiansen  effect,  which  is  particularly  pro¬ 
nounced  for  clouds  consisting  of  a  significant  number 
of  small  ice  crystals.  Based  on  a  line-by-line  equivalent 
solar  model,  we  show  from  spectral  curves  that  cloud 
reflection  increases  as  ice  crystal  sizes  become  smaller, 
but  the  cloud  absorption  increase  is  only  evident  for 
wavelengths  longer  than  about  2.7  /im.  The  ice  crystal 
shape  has  a  substantial  effect  on  the  cloud  reflection  and 
absorption  for  a  given  size;  more  complex  icc  particles 
reflect  more  solar  radiation.  Finally,  we  propose  a  con¬ 
trail  cirrus  cloud  model  consisting  of  a  combination  of 
bullet  rosettes  (60%),  hollow. columns  (30%),.andplates 
{’20%),  with  sizes  ranging  from  1  to  90  pm  in  association 
with  radiation  perturbation  studies. 


Introduction 

One  of  the  fundamental  objectives  of  Subsonic  Air¬ 
craft:  Contrail  and  Cloud  Effects  Spedd  Sta^  ^CUC- 
C£5S)  is  to  provide  the  necessary  and'sufl^ieHt  data 
to  better  determine  the  radiative  properties  of  contrail 
cirrus  and  their  potential  impact  on  climate.  Contrails 
tend  to  form  in  the  vikanity  of  the  tropopause  in  which 
the  ambient  temperature  allows  supersaturated  water 
vapor  emitted  from  jet  engines  to  effectively  convert 
into  ice  phase.  Resident  times  for  particles  may  vary 
from  days  below  the  tropopause  and  up  to  a  year  above 
the  tre^pause  in  the  lower  stratosphere.  Analysis  of 
air  traffic  patterns  illustrates  that,  m  general,  cruising 
altitudes  are  vdthin  contrail  cirrus  producing  levels.  A 
number  of  statistical  studied  have  indicated  a  significant 
correlation  between  jet  fuel  consumption  and  high  cloud 
frequencies  in  urban  airports  [Liov  et  of.,  1990;  Frankel 
et  at,,  1997].  In  particular,  the  data  from  Salt  Lake  City 
appear  to  support  the  connection  between  increasing 
jet  aircraft  operations  in  tfasjippeL  troppsphar^^ 
cloud  cover,  and  r^onal  climate  when  a  sharp  increase 
in  domestic  air  traffic  occurred  in  the  mid-l^s. 

Because  of  the  lack  of  in  stfu  micropbystcal  observa^ 
tions  in  contrails,  their  radiative  propWties  are  largely 
unknown.  The  unique  SUCCESS  field  experiment  car¬ 
ried  out  over  Kansas  in  April- May,  1995  provided  mi- 
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crophysical  measurements  on  the  size  and  shape  char¬ 
acteristics  ot  ice  crystals  that  were  not  available  in  the 
past.  Using  the  observed  ice  crystal  size  distributions 
and  shape  rectors  for  contrails,  we  have  perform^  scat¬ 
tering  and  absorption  calculations  based  on  a  unified 
theory  for  light  scattering  by  ice  crystals  introduced  in 
the  text.  Pertinent  results  for  the  extinction  coefficient 
smgle-scattering  albedo,  and  phase  function  for  contrail 
cirrus  are  presented.  Employing  a  lin©*by-line  equiva¬ 
lent  solar  model,  we  have  also  <fcterminea  the  spectral 
and  broadband  radiative  properties  of  contraul  cirrus. 


Ice  Crystal  Size  Distributions 

During  SUCCESS  on  May  4, 1996,  iceccystad  size  dis¬ 
tributions  for  the  contrails  were  measured  by  the  repli¬ 
cator  system  mounted  on  NASA^s  DC-8  that  tailed  a 
B^ing  767  bv  about  50  sec  with  aircraft  separation 
distance  of  11.5  km  over  northeast  Oklahoma.  This 
replicator  system  was  developc^JbyJh^  scientists  at  the 
Desert  Research  Institute  (DRI),  University  of  Nevada 
[Amott  et  aL,  1994J.  Two  representative  icc  crystal  size 
distributions,  associated  with  the  ambient  temperature 
and  dew  point  of  -fil.l’C  and  -62.9®C,  respectively, 
were  selected  for  the  present  study.  The  growth  of  ice 
crystals  in  the  contrails  took  about  50  sec.  The  icc  crys¬ 
tal  images  collected  by  the  DRI  replicator  system  show 
that  they  were  predominately  bullet  rosettes,  columns, 
and  plates  with  sizes  ranging  from  about  1  /<m  to  about 
90  um. 

In  addition  to  the  preceding  data,  we  have  also  ob¬ 
tained  two  ice  crystal  size  distributions  from  an  ex¬ 
periment  sponsored  by  the  Department  of  Energy  that 
was  carried  out  over  the  Cloud  and  Radiation  Testbed 
(CART)  site  located  in  Northern  Oklahoma  and  South¬ 
ern  Kansas  on  April  18  and  19,  1994.  On  18  April,  a 
contrail  was  studied  above  the  CART  site  when  the  Uni¬ 
versity  of  North  Dakota  Citation  aircraft  ra'penetrated 
its  own  contrail  after  6  mtn  of  growth  at  a  height  of 
13  km  and  a  temperature  of  -65.9®C.  As  shown  in 
Sassen  [19971,  the  DRI  replicator  indicated  minute,  sim¬ 
ple  plate  and  column  crystals  in  the  contrail.  The  data 
for  19  April,  on  the  other  band,  came  from  near  the  top 
(13.4  and  -69.4®C)  of  an  optically  thin  cirrostratus 
cloud  with  embedded  contrails,  which  contained  minute 
solid  ice  particles.  Since  this  cloud  generated  a  lunar 
corona  display,  as  contrails  frequently  do,  we  consider 
this  cloud  composition  as  a  proxy  oi  a  persisting  con¬ 
trail.  The  size  distribution  data  were  derived  from  the 
FSSP  device,  which  has  been  shown  to  yield  reliable 
data  when  appreciable  numbers  of  ice  crystals  larger 
than  50-100  arc  not  present  \Gayet  et  at.,  1996J. 

To  characterize  the  effect  of  ice  crystal  size  distriou- 
tions  on  the  radiative  transfer  results  and  to  account  for 
various  shapes  and  sizes,  we  define  the  mean  effective 
ice  crystal  size  io  terms  of  the  maximum  dimension,  D, 
in  the  form 

D,  =  j  Vn{D)dD/  j  An{D)dD  =  (1) 


Pipo- number  97OLO3}08.  whet«  V  and  A  at€  the  volume  and  projection  area  of  an 

OW4-,}3WW70U»50B$05.00  ice  crystal,  respectively,  fH^C  is  the  ice  water  content, 
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Figure  1.  Discretized  ice  crystal  size  distributions  for 
contrail  cirrus  50  sec  duration)  measured  by  the 
replicator  system  mounted  on  NASA's  DC-8  that  tailed 
a  Boeing  757  during  SUCCESS  on  May  4,  1996.  and  for 
a  contrail  and  a  cold  cirrus  6  min  duration)  measured 
by  FSSP  on  board  the  University  of  North  Dakota  Ci¬ 
tation  on  April  18  and  19,  1994  (upper  panel). 


Pi  is  the  bulk  ice  density,  and  At  is  the  projected  area 
per  unit  volume.  The  mean  effective  ice  crystal  size  so 
defined  is  now  directly  related  to  [WC^  a  prognostic 
parameter  in  GCMs.  Also,  by  using  the  volume  and 
projection  area,  the  irregular  ice  emtaJ  shape  can  be 
accounted  for.  Finally,  At  can  be  determined  directly 
from  the  ice  crystal  images. 

For  stngle*3cattering  calculations,  the  observed  ice 
crystal  size  distributions  are  discretized  in  eight  bins: 
1-5  3  ^m:  5-10  with  bin- 

center  atT:S'TZW7i0-20  ;im  with  bin-cciitcr  at  15  wm; 
20-30  pm  with  bin-center  at  25  ^m;  30-40  pm  with  bin* 
center  at  35  Mm;  50-70  Mm  with  bin-center  at  60  pm\ 
7CP90  Mm  with  biu-center  at  80  Mm,  Shown  in  Figure  1 
are  the  discretized  ice  crystal  size  distributions.  A  sub¬ 
stantial  number  of  small  ice  crystals  on  the  order  of  10 
Mm  is  seen  for  the  contrail  that  occurred  on  April  18, 
1994. 

Based  on  the  SUCCESS  replicator  data,  we  esti¬ 
mate  that  contrails  consist  approximately  of  50%  bullet 
rosettes,  30%  hollow  columns,  and  20%  plates.  Using 
these  combinations,  the  mean  effective  sixes  for  the  four 
ice  costal  size  distributions  are  4.9,  9.8,  15-9,  and  13.3 
Mm.  The  contrails  observed  on  April  18  and  19  contain 
smaller  ice  particles  than  those  ooserved  on  May  4. 


Phase  FuiiLtroii'and  Single-Scattering 
Properties 

Wc  have  used  the  Monte-Carlo/geometric-ray-tracina 
3U3Ptoach  to  determine  the  scattering,  absorption,  ana 
polarization  properties  of  plates,  solid  and  hollow  column: 
dendrites,  bullet  rosettes,  aggregates,  and  ice  crystals 
with  irre^lar  surfaces  whose  sizes  are  much  larger  than 
the  incident  wavelength  [TaAono  and  Liou,  1995].  For 
light  scattering  by  small  ice  crystals,  we  have  inno¬ 
vated  an  improved  geometric*ray*tracing  method  and 
a  finite-difference  time  domain  (PDTD)  method  [  Von^ 
and  Ltou,  1996a, b].  The  improved  model  is  applicable 
to  size  parameters  as  small  as  about  IS.  The  FDTD 
method  is  a  numerical  technique  for  the  solution  of  the 
Maxwell  equations  in  time  and  space  using  appropriate 
absorbing  boundary  conditions  It  is  considered  to  be 


Figure  2*  Phase  functions  for  0.7  and  3  7  Mm  wave¬ 
lengths  using  a  contrail  cirrus  model  consisting  of  50% 
bullet  rosettes,  30%  hollow  columns,  and  20%  plates 
(left  panel).  The  effect  of  ice  crystal  surface  rou^ness 
on  phase  function  is  shown  in  the  right  panel  using  bul¬ 
let  rosette  aa  an  example 


the  exact  numerical  solution  for  light  scattering  by  par¬ 
ticles,  as  verified  by  the  exact  theoretical  solutions  for 
long  circular  cylinders  and  spheres.  Because  of  numer¬ 
ical  round-off  errors  and  required  computer  CPU  time, 
the  FDTD  method  can  only  be  applied  to  ice  crystal 
size  parameters  smaller  than  about  10-15.  By  unify¬ 
ing  tne  modified  geometric  optics  and -ED'FD. methods^ 
for  large  and  smaQ  size  parameters,  respectively^  we  arc 
now  in  a  position  to  determine  the  basic  «ettenng,  ab¬ 
sorption,  and  polarization  properties  of  ice  cry^als  of 
any  sizes  and  snapes  that  can  be  defined  numerically. 

The  sintfle-flcattering  parameters  in  terms  of  phase 
function,  single-scattering  albedo,  and  extinction  cross 
section  are  computed  for  the  solar  spectrum  covering 
0.2  to  5  pm  for  about  200  wavelengths.  The  diagram 
in  the  left  panel  in  Figure  2  show  the  phase  functions  tor 
the  wavelengths  of  0  7  and  3  7  pm^  Substantial  differ¬ 
ences  in  the  phase  function  for  the  four  cases  are  noted 
at  3.7  pm  at  which  ice  is  strongly  absorptive  and  the 
ice  crystal  size  has  a  significant  impact  on  the  singly 
scattering  feature.  For  the  case  of  4/18/94,  the  halo 
peaks  and  backscattering  for  the  0.7  pm  reduce  sub¬ 
stantially  because  of  smaul  ice  crystal  sizes. 


Figure  3.  Asymmetry  factor,  single-scattering  albedo, 
and  extinction  coefficient  as  functions  of  wavelength 
from  0.2  to  6  fim.  The  minima  located  at  2.85  ftm 
are  the  well-known  Christiansen  effect. 
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From  the  replicator  data,  it  appears  that  some  ice 
crystals  may  not  be  exactly  hexa^nal  in  shape  and  that 
some  may  have  rough  surfaces.  FollowinE  Takano  and 
Liou  (1995],  we  have  used  a  Monte-Cano/geometric- 
ray^tracing  approach  to  account  for  the  Frcroclian  inter¬ 
action!  between  rays  and  rough  surfaces  of  an  ice  crystal 
defined  by  the  fitet-order  Gram  Charlier  Series  distribu¬ 
tion  [Cox  and  AfunJb,  19541.  The  diagrams  in  the  right 
panel  of  Fig.  2  illustrate  the  phase  function  differences 
Between  smooth  and  rough  surfaces  using  bullet  rosettes 
as  an  example.  The  smooth-surface  bullet  rosettes  pro¬ 
duce  halo  peaks  at  22^.  46^,  and  10^,  as  well  as  the 
backscattenng  peak  ana  a  maximum  feature  at  about 
160^  at  0.7  am.  All  of  these  peak  features  are  smoothed 
out  when  the  roughness  on  the  surface  is  incorporated 
in  the  light  scattering  calculations.  For  strongly  absorp- 
tive  cases,  the  effect  of  surface  roughness  on  the  phase 
function  for  randomly  oriented  ice  crystals  is  small  be¬ 
cause  the  scatt^ing  teatures  are  generated  primarily  by 
external  reflections. 

Figure  3  shows  the  asymmet^  factor,  single-scattering 
albedo,  and  extinction  coefficient  based  on  the  shape 
model  of  50%  bullet  rosettes,  30%  hollow  columns,  and 
20%  plates.  The  extinction  coefficients  show  little  vari¬ 
ation,  except  for  a  minimum  in  the  2.85  am  region, 
the  so-callM  Christiansen  effect.  This  effect  occurs 
when  the  real  part  of  the  refractive  index  approaches 
1  while  the  corrMponding  imaginary  part  is  sul^tan- 
tially  larger,  leading  to  the  domination  of  absorption  in 
light  attenuation.  It  is  particularly  pronounced  when 
ice  particles  are  small,  such  as  the  4/18/94  caw.  The 
sin^e-scattering  alMo  ^s6  displays  a  strong  minimum 
in  %e  2.85  /<m  region  with  values  much  less  than  0.5. 
When  absorption  is  substantial,  the  scattered  energy 
is  primarily  contributed  by  djjBEraction  in  the  forwatd 
directions.  For  this  reason,  maximum  values  of  the 
asymmetry  factor  are  noted  around  3  /im*  Introduc¬ 
ing  the  roughness  of  ice  crystal  surfaces  does  not  eig- 
^cantly  meet*  the  extin^on  coefficient  and  single- 
scattering  albedo  patterns.  It  should  also  be  noted  that 
eoot  particles  can  be  activated  as  ice  nuclei  in  contrail 
formation  and  the  effect  on  the  scattering  and  absorp¬ 
tion  of  ice  particles  is  dependent  on  their  relative  sises. 


Spectral  and  Broadband  Solar  Albedos 

We  have  developed  a  new  solar  radiation  rnodel  with 
a  high  epcctral  resolution  similar  to  Une-by-Une  calcu¬ 
lations  in  the  thermal  IR,  hereafter  referred  to  as  the 
line-by-line  equivalent  (LBLE)  solar  model  The  spec¬ 
tral  resolution  has  two  options:  1  with  10 

and  50  with  30  i?*s,  where  g  denotes  the  cumula¬ 
tive  probability  function  in  the  corrdated  k-distribution 
metb^.  The  correlated  k^coeffidents  for  these  spectral 
intervals  for  H5O  from  2,000  to  21,000  fO.5-5  /im) 

are  determined  from  the  updated  1996  HITRAN  data 
baaed  on  the  method  developed  by  Fu  and  Liou  [1993]. 
The  correlated  k-ooefficients  for  the  2.0  and  2  J  /imCOs 
bands  are  also  derived,  in  which  overlaps  between  H^O 
and  COs  lines  arc  accounted  for  by  means  of  the  mul¬ 
tiplication  rule.  Absorption  due  to  O3  and  Og  bands 
and  Rayleigh  scattering  contributions  are  accounted  for 
baaed  on  the  conventional  method.  In  addition  to  cloud 
particles,  we  have  also  compiled  the  single-scattering 
properties  of  six  typical  aerosol  in  connection  with  the 
lBLE  solar  model. 

The  **exact^  adding/doubling  radiative  transfer  pro- 
mm  including  full  otokes  parameters  developed  by 
Tdkano  and  Liou  [1989]  is  used  to  calculate  the  transfer 
of  monochromatic  solar  radiation  in  a  vertically  inho- 
mogeneou!  atmosphere  which  is  divided  into  61  lay^. 
For  wavelengths  betwwn  3.7  to  5  ^m,  thermal  emission 


Figure  4.  Spectral  reflection  and  absorption  for  con¬ 
trail  cirrus  as  a  function  of  the  mean  effective  ice  crystal 
size  for  two  solar  zenith  angles.  The  contrail  cirrus  were 
located  at  12  km  over  a  land  surface  with  a  Lambertian 
albedo  of  0.2  and  with  an  IW P  of  1  gm“^. 


contributions  to  the  solar  flux  transfer  are  accounted  for 
by  adding  the  thermal  emission  part,  (I  -  w)xBx(r), 
to  the  adding/doubling  method,  where  u;  is  the  single- 
scattering  albedo  and  ;rFx(T)  is  the  Planck  flux  for  a 
given  temperature.  The  thermal  part,  although  small, 
has  not  been  accounted  fori^orfaowfaaad  solar  flux  cal¬ 
culations  and  could  be  a  noticeable  en^gy  source  in 
the  upper  part  of  the  atmosphere;  a  subj^  for  future 
exploration.  The  spectral  solar  radiative  transfer  pro¬ 
gram  also  includes  the  options  of  using  the  detailed  so* 
lar  irrs^iance  data  averaged  over  appropriate  spectral 
resolutions. 

Piwrc  4  illustrate  the  spectrgJ  and  edjsorp- 

tion  for  the  contrail  layer  .sSLja..huiction  of  the  mean 
effective  ice  costal  size  for  two  solar  zenith  angles  of 
30*  and  75®.  Tb  economize  the  computations,  we  have 
used  the  spectral  version  of  SO  containing  30 
From  observations,  the  contrails  were  located  at  about 
12  km  over  a  laind  surface  which  is  assumed  to  have  a 
Lambertian  albedo  of  0.2.  In  the  calculations,  we  have 
employed  an  ice  water  path  (IWF)  of  1  cm”*.  The 
visible  optical  depth  for  an  ice  crystal  cloud  can  be  ex- 
preu^  in  terms  of  both  IWP  and  ice  crystal  size  in 
the  form  [JFb  and  Liou,  1993] 

T,  =  I{VP{c^b/Pe).  (2) 

where  a  and  6  are  certain  coefficients.  Thus,  for  a  given 
IWP,  the  cloud  r^ection  is  inversely  proportional  to 
the  mean  effective  ice  cryst^  size  P^i  as  illustrated  in 
Fig.  4.  THe  wafer" vapor  abR&ption  oahds  located  at 
3.2,  2.7,  1,87,  1,38,  l.l,  0.94>-e;82,  and  0.72  are 
clearly  shown,  as  is  the  4.3  am  carbon  dioxide  band. 
The  large  reflection  in  the  UV  region  b  primsuily  pro¬ 
duced  by  Rayleigh  scattering.  For  the  solar  zenith  angle 
of  75®,  substantial  reflection  Is  evident  for  wavelengths 
shorter  than  about  2.5  fim,  particularly  for  smaller 
The  effect  of  ice  crystal  size  on  cloud  absorption  appears 
only  for  wavelengths  greater  than  2.8  at  whidi  ab¬ 
sorption  due  to  ice  particles  becomes  more  significant 
than  absorption  due  to  water  vapor  and  carbon  dioxide. 

The  broadband  reflection  and  absorption  for  contrail 
layers  as  a  function  of  IWP  arc  shown  in  Figure  5,  The 
left  and  right  panels  illustrate  the  effects  of  ice  crystal 
size  and  shape  on  the  cloud  radiative  prwerties,  respro- 
tivcly.  In  the  calculations,  we  have  used  a  solar  zenith 
angle  of  30®  and  a  surface  albedo  of  0.2.  Other  input 
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Figure  5.  Broadbuid  reflection  and  a^orption  for 
contrail  cirrus  as  a  function  of  IWP  to  illustrate  the 
effects  of  ice  crystal  rise  (left  panel)  and  shape  (right 
panel).  The  solar  zenith  angle  used  is  30". 


parameters  are  the  same.  Oloud  reflection  for  a  given 
tWP  is  clearly  dependent  on  ice  crystal  size;  smaller 
ice  crystals  reflect  more  solar  radiation  than  larger  ice 
crystals.  When  IWP  is  smaller  than  ateut  20  jm”  , 
cloud  absorption  also  reduces  when  the  ice  crystal  size 
increases.  For  tWP  larger  than  about  20  jm  .  cloud 
absorption  for  the  esse  of  4.9  /im  bewimes  the  sraailCTt 
because  it  is  now  determined  by  the  siMle^scattcring 
albedo  rather  than  the  optied  depth.  To  investigate 
the  shape  effect,  we  have  used  a  Dt  of  15.9  #im.  The 
effect  Of  rou^ness  of  bullet  rosette  sutfwes  cm  cloud  re¬ 
flection  and  absorption  is  insignificant.  Cloud  reflection 
values  become  progresavely  smaller  for  hollow  columns, 
plates,  and  equal-area  spheres  because  their  asymme¬ 
try  factors  become  increasingly  larger  to  allow  str^e* 
forward  scattering  to  occur.  Cloud  absorption  vstues 
for  plates,  hollow  columns,  and  equal-area  sphew  be¬ 
come  larger  bwause  of  the  progressively  smaller  single- 
scattering  albedo  values. 


ice  crystal  shape  has  a  substantial  effect  on  the  reflec¬ 
tion  and  absorption  calculations.  For  a  given  ice  crystal 
sise,  cloud  reflection  values  decrease  progressively  wheu 
bullet  rosettes,  hollow  columns,  plates,  and  equal-area 
spheres  are  employed  in  the  calculations. 

In  order  to  investigate  reliably  the  radiative  and  cli¬ 
matic  effects  of  contrail  cirrus,  an  appropriate  ice  crys¬ 
tal  model  in  terms  of  size  «id  shape  is  requited.  We 
have  proposed  a  combination  of  bullet  rosettes,  hollow 
colmnns^  and  plates  with  sises  tangms  from  I  to  90 
iim  for  contrail  cirrus  in  association  with  radiation  per¬ 
turbation  studies.  Finally,  the  subject  of  the  finite  and 
inhomogeneous  nature  of  contrail  cirrus  on  the  radiative 
effects  is  one  that  requires  further  investigation. 
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Summary 

We  have  analysed  two  ice  crystal  size  distributions 
for  contrail  cirrus  collected  dunng  SUCCESS  and  ob¬ 
tained  two  addUtional  distributions  from  an 
that  was  carried  out  earlier  over  the  DOE  CART  site. 
Using  these  four  observed  ice  costal  sire  distributions, 
we  have  performed  the  scattering  and  absorption  cal¬ 
culations  to  determine  the  single-scattering  and  pha^ 
function  properties  of  ice  particles  that  were  present  in 
coutraiU  for  wavelengths  from  0.2  to  5  Based  on 
the  SUCCESS  in  situ  observed  data,  we  prop^^tpee 
crystal  shape  model  for  contrails  consisting  of  50%  bul¬ 
let  rosettes,  30%  hollow  columns,  and  20%  plates  in 
conjunction  with  radiation  calculations. 

For  strongly  absorptive  cases,  the  ice  crystal  sur¬ 
face  roughness  is  not  a  significant  factor  in  determin¬ 
ing  the  phase  function  pattern.  The  extinction  coef- 
fiaent  Shows  a  minimum  feature  at  2.85  pm  referred 
to  as  the  Christiansen  effect,  pattlcutoly  pronounce 
when  small  ice  crystals  are  present.  The  correspond¬ 
ing  sinde-scattef  ing  albedo  in  this  spectral  region  shows 
v^ues  much  less  than  0.5.  Increasing  the  ice  crystal  size 
increase  the  asymmetry  factor,  whereas  the  reverse  is 
true  for  the  single-scattering  albedo. .  ,  ,  . 

Based  on  a  LBLE  solar  model,  the  dependence  of 
cloud  reflection  on  the  mean  effective  ice  crystal  size  for 
a  given  IWF  is  clearly  illustrated  in  the  entire  solar 
9p^trum.  For  a  given  IWP^  smaller  ice  crystals  re- 
jfl^  more  solar  radiation  than  larger  ice  crystals.  The 


154-169,  Springer  Verlag,  1990. 

Sasten,  K.,  Contrail-cirrus  and  their  potential  for  regional 
climate  change,  BuiL  Amer,  Meteor,  Soc.,  7$,  1885-1903, 
1997.  ^  . 

Takano,  Y-,  and  K-  N.  Liou,  Solar  radiative  transfer  m  cirrus 
clouds.  II:  Theory  and  computation  of  multiple  scattering 
In  an  anisotropic  medium,  J.  Atmos,  Sci,,  46, 20-36,  1989, 
Takano,  Y.,  and  K.  N  Liou,  Radiative  transfer  in  rirrus 
clouds.  Ill;  Li^t  scatttaing  by  irregular  ice  crystals,  J. 
Atmos,  Sci,,  818-837,  1995- 
Yang,  P.,  and  K.  N-  Liou,  Finite-difference  time  domain 
. ..  method  for  li^t  scattering  by  small  ice  crystals  in  three- 
dimensional  space,  /.  Opt.  Soc,  Am.  A,,  IS,  2072-2085, 
1996a. 

Yang,  P-1  and  K.  N.  Liou,  Georaetric-opticsdntcgral- 
equation  method  for  light  scattering  by  nonspherical  ice 
crystsJb,  Appt  Opt,  SS,  6568-6584,  1996b. 

K-  N.  Liou,  P.  Yang,  and  Y.  Tbkono,  Department  of 
Atmospheric  Sciences,  UCLA,  Los  Angeles, CA  90095-1565 
K.  Sassen,  Department  of  Meteorology,  University  of 
Utah,  Salt  Lake  City,  UT  84112  .  .  .  .  , 

T.  P.  Charlock,  Atmospheric  Science*  Diviaion,  Langley 
Researdi  Ceatcr/NASA,  Hampton,  VA  23665 
W.  P.  Amott,  Atmospheric  Science*  Center,  Desert 
Researdi  Institute,  BetX),  NV  89506 

(Received  July  10,  1997;  revised  November  20,  X997: 
accepted  Novembw  25,  1997.) 


AIRBORNE  RETRIEVAL  OF  CIRRUS  CLOUD  OPTICAL  AND 
MICROPHYSICAL  PROPERTIES  USING  ARES 
5. 1-5.3  AND  3.7  |im  CHANNEL  DATA 


S.  C.  Ou,  K.  N.  Liou,  P.  Yang,  and  P.  Rolland 
Department  of  Atmospheric  Sciences,  UCLA 
Los  Angeles,  California 

T.  R.  CaudiU 

Air  Force  Research  LaboratoryA/^SBE 
Hanscom  AFB,  Massachusetts 

J.  Lisowski 

SCITEC,  Princeton,  New  Jersey 
B.  Morrison 

Aeromet,  Inc.,  Tulsa,  Oklcihoma 


Abstract 


We  present  an  airborne  retrieval  algorithm  to  infer  cirrus  cloud  temperature,  optical  depth  and 
mean  effective  sizes  using  the  Airborne  Remote  Earth  Sensing  System  (ARES)  5. 1-5.3  and  3.7  pm 
channels  data.  The  algorithm  development  and  the  selection  of  the  channels  are  based  on  the  principle 
and  parameterization  of  radiative  transfer  involving  cirrus  clouds  and  the  associated  atmospheric  and 
surface  properties.  It  has  been  applied  to  the  ARES  data  collected  over  the  western  Boston  area  on 
September  16,  1995.  Validation  of  the  retrieved  parameters  was  carried  out  using  the  collocated  and 
coincident  ground-based  8.6-mm  radar  data  and  ice  crystal  size  distribution  measurements  obtained  from 
the  2D-P  probe  on  board  the  High  Altitude  Reconnaissance  Platform  (HARP).  We  show  that  the 
retrieved  cirrus  cloud  temperature,  mean  effective  ice  crystal  size,  and  optical  depth  match  closely  with 
those  derived  from  the  observations. 
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1 .  Introduction 

Established  cirrus  cloud  retrieval  methods  based  on  imaging  radiometric  data  generally  use  the 
visible,  window  (e.g.,  3.7,  10.9,  and  12.0  urn),  and  water  vapor  absorption  (6.3  pm)  channel  radiances 
that  are  available  from  either  operational  satellites  or  research  aircraft  (Liou  et  al.  1990;  Minnis  et  al. 
1990,  1993;  Baum  et  al.  1994;  Ou  et  al.  1993,  1995).  Remote  sensing  of  the  cirrus  cloud  parameters 
using  radiances  from  other  spectral  channels  has  also  been  recently  explored.  Wielicki  et  al.  (1993) 
investigated  the  possibility  of  using  the  0.83  and  1 .65  pm  radiances  from  Eandsat  to  determine  ice  crystal 
particle  sizes.  Gao  and  Kaufman  (1995)  reported  a  new  technique  for  detecting  thin  cirrus  clouds  using 
the  1.38  and  1.88  pm  water  vapor  absorption  channels  radiances  from  the  Airborne  Visible  Infrared 
Imaging  Spectrometer  (AVIRIS),  based  on  the  fact  that  the  low-level  water  vapor  absorption  provides  a 
dark  background  for  the  high-altitude  cirrus  reflection  to  be  detectable.  The  1.6  and  2.2  pm  channel 
radiances  from  the  Moderate  Resolution  Imaging  Spectrometer  (MODIS)  and  the  MODIS  Airborne 
Simulator  (MAS)  are  also  potentially  useful  channels  for  future  developments  of  cirrus  cloud  retrieval 
algorithms  (King  et  al.  1992,  1995;  Tsay  et  al.  1996),  because  both  wavelengths  are  sensitive  to 
variations  in  ice  crystal  size  distributions.  Moreover,  the  selection  of  appropriate  channels  for  reliable 
cirrus  cloud  detection  and  retrieval  will  be  a  major  consideration  in  the  development  of  the  National 
Polar-orbiting  Operational  Environmental  Satellite  System  (NPOESS)  project. 

In  this  paper,  we  present  an  airborne  retrieval  algorithm  to  infer  cirrus  cloud  temperature,  optical 
depth,  and  ice  crystal  mean  effective  size  using  the  ARES  IR  spectrometer  data.  The  retrieval  method 
utilizes  for  the  first  time  the  5. 1-5.3  pm  water  vapor  band  to  determine  cirrus  cloud  temperatures  and  IR 
emissivities.  from  which  the  associated  optical  depths  can  be  determined  from  radiative  transfer 
parameterizations.  In  the  scheme,  the  ice  cry  stal  mean  effective  size  is  derived  by  matching  the  sum  of 
the  computed  3.7  pm  solar  reflected  and  IR  emitted  radiances  with  the  ARES  Ch.  42  (3.755.  pm) 
radiances.  The  selection  of  these  retrieval  channels  is  based  on  the  examination  of  the  computed 
weighting  functions  of  a  typical  midlatitude  summer  atmosphere  along  with  the  inspection  of  ARES 
channel  radiance  images.  We  have  applied  the  retrieval  algorithm  to  the  ARES  data  collected  during  a 
mission  flight,  which  was  carried  out  on  board  the  WB-57F  aircraft  on  September  16,  1995  over  the 
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western  Boston  area.  The  retrieved  cloud  parameters  are  compared  with  the  in-situ  2D-probe 
measurements  obtained  from  the  Aeromet's  HARP,  as  well  as  the  cloud  heights  derived  from  the  ground- 
based  8.6-mm  radar  returns. 

Cloudy  radiances  have  been  measured  by  the  IR  imaging  spectrometer/radiometer  on  board  the 
ARES  Programs'  WB-57F  high-altitude  research  aircraft  operated  by  the  NASA  Johnson  Space  Right 
Center.  The  WB-57F  is  capable  of  reaching  altitudes  in  excess  of  60,000  ft,  and  carries  a  crew  of  two, 
permitting  instrument  operation  by  a  dedicated  flight  crew.  The  sensor  was  mounted  on  the  lower 
fuselage  surface,  with  downward  viewing  while  in  level  flight.  The  ARES  sensor  for  the  present  study 
has  been  operated  as  a  75-channel  imaging  spectrometer.  It  is  operated  in  a  pushbroom  mode  with  45- 
pixel  array  pointing  in  the  nadir  direction  taking  data  at  the  standard  rates  of  10.  20,  40,  or  80  scans/sec. 
The  optical  elements  of  .ARES  provide  a  pixel  EFOV  (instantaneous  field  of  view)  of  slightly  more  than  1 
mrad,  which  corresponds  to  about  15  m  ground  resolution  for  a  flight  altitude  of  15  km.  The  spacing  of 
the  ARES  spectral  wavelengths  appears  to  depend  on  the  optical  arrangement  with  the  channel  band 
widths  that  vary  between  0.02  and  0.07  pm.  Further  details  of  the  ARES  can  be  obtained  from  the 
internet  web  site  :  http://cahill.plh.af.mil;8600/PDAC.html. 

For  validation  of  the  retrieved  cloud  temperature/height,  we  have  obtained  the  cloud  boundary 
height  data  from  the  collocated  and  coincident  ground-based  8.6-mm  TPQ-1 1  radar  measurements.  For 
comparison  of  the  retrieved  cloud  optical  depth  and  mean  effective  sizes  with  in-situ  measurements,  we 
have  also  acquired  the  meteorological  and  ice  microphysical  data  collected  from  the  HARP,  which  is  a 
Leaijet  36A,  flying  in  conjunction  with  the  ground-based  radar  and  the  ARES  measurements.  The  HARP 
carried  instruments  for  taking  both  in-situ  and  remote  atmospheric  measurements.  The  in-situ  instruments 
include  Forward  Scattering  Spectrometer  Probe  (FSSP),  ID-C  and  2D-C  probes  for  particle  size 
distribution  measurements.  Ambient  temperatures  were  measured  by  the  Rosemount  Total  Temperature 
Sensor  with  an  uncertainty  of  ±1  K. 

This  paper  is  organized  as  follows.  In  Section  2,  we  describe  the  development  of  the  remote 
sensing  algorithm  for  the  retrieval  of  cirrus  cloud  temperature  (height),  optical  depth,  and  ice  crystal  size 
principally  using  the  ARES  5. 1-5.3  and  3.7  pm  channels.  In  Section  3,  results  of  the  application  of  this 
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relneval  scheme  to  the  ARES  data  collected  during  September  16,  1995,  are  presented  along  with  the 

validation  efforts  using  the  ground-based  radar  and  airborne  in-sin.  measurements.  Finally,  a  summary  Is 
given  in  Section  4. 

2.  Principle  of  Retrieving  Cirrus  Parameters  Using  ARES  Thermal  IR  Channels 

The  present  retrieval  method  for  inferring  cirrus  temperature  and  emissivity  is  based  on  the 
numerical  solution  of  the  following  set  of  algebraic  equations  (Liou  etal.  1990;  Ou  etal.  1993,  1995)  • 

Ri  =  ( 1  -  £( )  Rji  +  Ei  B|  (Tc),  ^  j  ^ 

Ei  =  1  -  exp(-kiT),  ^2) 

where  R,  is  the  upwelling  radiances  at  the  top  of  the  cirrus  cloudy  atmosphere  for  the  ith  selected  ARES 
channel;  R^j  is  the  corresponding  radiances  reaching  the  cloud  base  for  the  ith  channel;  Ej  is  the  cloud  IR 
emissivities,  T^  is  the  mean  cirrus  temperature;  (Te)  is  the  Planck  intensities  at  Tg,  x  is  the  visible 
optical  depth  and  ki  is  the  effective  extinction  coefficient  that  has  been  determined  from  the  calculated 
single-scattenng  properties  based  on  a  number  of  measured  ice  cry  stal  size  distributions.  The  behavior  of 

the  cloud  emissivity  in  terms  of  its  dependence  on  wavelength  and  cloud  parameters  has  been  examined 
previously  (Liou  et  al.  1990). 

By  straightforward  algebraic  analyses,  Eqs.  (1)  and  (2)  for  a  pair  of  channels  (i=l  and  2)  can  be 
combined  into  a  single  equation  in  the  form 

[  R2  -  BzCTc)  ]  /  [  Ra2  -  B,(Tc)  ]  =  {  [  R,  -  B,(Tc)  ]  /  [R„  -  B,(Tc)  ]  } .  (3) 

If  Ra,  and  R,,  are  known.  B:(Tc)  and  B,(Te)  are  approximately  linearly  correlated,  and  k:/k,  ~  1.0,  then 

Eq.  (3)  becomes  a  linear  equation,  whose  numerical  solu'ion  can  be  obtained  in  a  relatively 
straightforward  manner. 

The  effeerive  excinction  coeff, dents  for  the  retrieval  channels  have  been  computed  using  the 
adding/doubling  program  in  conjunction  with  the  development  cf  this  remote  sensing  algorithm.  Eleven 
ice  crystal  size  distributions  from  aircraft  in-sim  measurements  have  been  chosen  (Fu  and  Liou  1993). 
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Their  mean  effective  ice  ctystai  sizes  range  between  23.9  and  123.6  Min.  The  stngle-scattering  propenies 
inciuding  phase  functions,  singie  scattering  albedos,  and  extinction  coefficients  have  been  computed  for 
each  size  distribution  and  for  each  of  the  selected  ARES  channels  using  the  unified  theoiy  for  light 
scattering  by  ice  crystals  (Yang  and  Liou  1996a,  b).  To  optimize  the  computational  effort,  the  forward¬ 
scattering  peak  of  each  phase  function  was  tnincated  following  the  pmcedure  developed  by  Takano  and 
Lion  (1989).  Radiances  for  each  channel  and  for  a  range  of  optical  depths  and  ice  crystal  sizes  were 
computed  using  the  aforementioned  adding/doubling  program.  The  cirrus  cloud  eraissivity,  which  is 
defined  as  the  ratio  of  the  upwelling  radiance  to  the  cloud  Planck  function,  was  obtained  as  functions  of 
the  optical  depth,  mean  effective  size,  and  channel  wavelengths.  Finally,  the  effective  extinction 
coefficients  for  each  channel  and  for  each  size  distribution  was  obtained  as  the  slope  of  the  linear  curve 
fitting  of  -ln(  l-a  vs,  t.  Figure  1  shows  the  effective  extinction  coefficient  (k-value)  as  a  function  of  the 
mean  effective  size  for  eight  wavelengths.  For  all  wavelengths,  the  k-values  increase  with  the  mean 
effective  size.  This  functional  dependence  of  the  k-values  has  been  incorporated  into  the  retrieval 
program.  More  importantly,  the  diagram  indicates  that  the  values  of  the  three  pairs  of  k-ratios  -  k(5.l89 
Mm)/k(5.1 16  pm).  k(5.261  Mm)/k(5.l  16  pm),  and  k(5.332  pm)ft(5.l  16  pm)-depend  weakly  on  the  mean 
effective  size  and  are  all  close  to  1.0  for  the  range  of  the  mean  effective  size  considered  here.  Thus,  they 
are  fixed  at  1 .0  in  the  numerical  solution  of  Eq.  (3). 

In  order  to  infer  the  cirrus  emissivity,  e,  and  temperature,  T.,  from  multiple  radiance 
measurements  in  cloudy  conditions,  the  upwelling  radiances  at  the  cloud  base.  R„,  must  be  known.  In 
principle,  it  can  be  generated  from  forward  radiation  transfer  calculations  if  the  surface  temperature, 
atmospheric  temperature  and  humidity  profiles,  and  cloud  base  height  are  given.  However,  a  more  direct 
and  reliable  method  to  determine  R„  would  be  to  find  cleat  column  radiances  from  the  available  data. 

Assuming  that  the  water  vapor  absorpiion/emission  effects  are  negligible  above  cions  clouds,  we 
may  approximate  the  cloud-base  upwelling  radiance  using  the  measured  clear  radiance  at  the  ARES  flight 
level.  This  approximation  has  been  applied  to  the  retrieval  of  clous  cloud  parameters  using  the  AVHRR 
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3.7  and  10.9  nm  window  bands  (Ou  et  al.  1993)  and  to  the  airborne  cirrus  retrieval  using  the  6.5  and  10.5 
Hm  radiometer  data  (Liou  et  al.  1990).  Once  is  determined,  Eq.  (3)  can  be  solved  by  a  numerical 
iteration  method.  Since  there  are  two  independent  parameters,  N-1  pairs  of  the  N  (>  1)  channel  data  are 
used  in  the  retrieval  program.  N-1  sets  of  the  retrieved  temperature  and  emissivity  are  then  averaged  to 
obtain  their  respective  mean  values. 

Upon  retrieval  of  the  IR  emissivity  and  the  cloud  temperature  using  the  thermal  IR  channels,  we 
can  obtain  the  mean  effective  size  using  the  observed  3.7  pm  channel  radiances.  We  first  assume  a  trial 
ice  crystal  size,  from  which  the  visible  optical  depth  can  be  obtained  from  Eq.  (2).  The  3.7  pm  thermal 
IR  radiance  can  be  calculated  according  to  Eq.  (1).  The  cloud-base  upwelling  3.7  pm  IR  radiance  can  be 
obtained  by  converting  the  4.75  pm  dear-pixel  radiance  into  its  equivalent  Planck  function  at  3.7  pm, 
since  both  channels  are  in  the  window  spectral  interval.  The  3.7  pm  solar  reflectance  is  determined  from 
look-up  tables,  which  contains  the  computed  value  of  the  3.7  pm  solar  reflectance  as  functions  of  the 
visible  optical  depth  and  the  mean  effective  size.  These  tables  are  generated  for  the  same  six  size 
distributions  as  described  in  Rao  et  al.  (1995)  and  for  appropriate  solar-viewing  geometries,  using  the 
adding/doubling  method  developed  by  Takano  and  Liou  (1989).  Finally,  the  total  3.7  pm  radiance  is 
calculated  as  the  sum  of  the  thermal  IR  and  solar  radiances.  The  radiance  values  thus  obtained  are 
subsequently  compared  with  the  observed  3.7  pm  radiances.  The  mean  effective  ice  crystal  size  can  then 
be  obtained  iteratively  by  adjusting  the  trial  value  based  on  the  difference  between  the  computed  and 
observed  radiances.  Pixels  having  a  retrieved  mean  effective  size  greater  than  123.6  pm  or  less  than  23.9 
pm  are  discarded.  Figure  2  gives  a  schematic  flow  diagram  for  the  retrieval  of  cirrus  cloud  temperature, 
emissivity  (optical  depth)  and  mean  ice  crystal  size. 

3.  Application  of  the  Retrieval  Scheme  to  the  ARES  Spectrometer  Data 

The  preceding  cirrus  cloud  retrieval  scheme  is  applied  to  the  ARES  spectrometer  data  obtained 
during  the  WB-57F  flight  that  was  conducted  on  September  16,  1995.  The  flight  has  been  divided  into 
eight  tracks.  Data  collected  for  Tracks  1.  3,  5,  and  8  were  analyzed.  These  four  tracks  are  straight,  and 
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spectrometer  data  were  recorded  in  nadir-pointing  pushbroom  mode.  The  flight  time  for  each  track  is: 
2055-2100  UTC  (Track  1).  2135-2140  UTC  (Track  3),  2212-2217  UTC  (Track  5),  and  2308-2313  UTC 
(Track  8).  The  cruise  portion  of  each  track  was  at  an  altitude  of  15  km  along  the  southeast-northwest 
direction,  crossing  over  the  Hanscom  Air  Force  Base  at  Bedford,  Mass.,  where  the  ground-based  8.6-mm 
radar  was  located.  During  the  period  of  Tracks  1  and  3,  HARP  was  also  flying  at  the  cirrus  cloud  level 
below  the  WB-57F.  Between  2030  and  2040  UTC  (case  1),  and  between  2147  and  2157  UTC  (case  2), 
HARP  flew  in  a  race  track  pattern  which  was  no  more  than  30  km  from  the  Hanscom  Air  Force  Base.  Its 
flight  tracks  are  approximately  collocated  and  coincident  with  the  route  of  the  ARES  Tracks  1  and  3, 
respectively.  For  both  cases,  the  flight  altitude  was  at  about  8  km,  close  to  the  top  of  the  cirrus  cloud 
vertical  domain. 

3. 1  ARES  Channel  Spectral  Properties 

Implementation  of  the  preceding  retrieval  scheme  requires  understanding  of  the  ARES  channel 
spectral  properties  and  the  appropriate  selection  of  useful  ARES  channel  data.  Figure  3  presents  the 
upwelling  radiances  associated  with  solar  reflection  and  thermal  IR  emission  in  the  2-6.4  pm  region  to 
assess  the  relative  importance  of  the  two  radiation  sources.  For  the  solar  reflection,  we  have  used  the 
updated  solar  irradiance  data  provided  for  our  use  by  G.  Anderson  at  the  Air  Force  Geophysics 
Directorate  (Anderson  1995).  The  results  are  shown  for  the  cosine  of  two  solar  zenith  angles  and  for  a 
surface  reflectance  of  0.1.  For  the  thermal  IR  emission,  Planck  function  curves  for  three  temperatures 
covering  the  range  of  relevant  earth's  atmospheres  are  shown.  It  is  clear  that  for  X,  <  3  pm,  the  solar 
reflection  is  the  primary  contributor  to  radiation,  while  for  X  >  4.5  pm,  the  thermal  IR  emission  is  the 
predominant  radiation  source.  For  3  pm  <  X  <  4.5  pm,  the  relative  importance  of  the  solar  reflection  and 
thermal  IR  emission  depends  on  the  sun's  position,  the  surface  albedo,  and  the  atmospheric 
thermodynamic  state.  For  the  radiative  transfer  simulation  and  the  remote  sensing  algorithm 
development,  both  components  must  be  accounted  for  in  this  intermediate  spectral  range. 
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Based  on  the  theory  of  atmospheric  gaseous  absorption,  several  absorption  bands  exist  in  the  2- 
6.4  pm  spectral  interval  (Goody  and  Yung  1989).  These  include  the  2.7  pm  H2O  band  produced  by  v, 
and  V:  fundamentals,  the  3.2  pm  H2O  band  associated  with  Zvi,  the  4.3  pm  CO2  band  due  to  the  v., 
fundamental  and  P  and  R  branches,  and  the  HjO  6.3  pm  vibrational  (V2)-rotational  band.  There  are  three 
windows  located  at  the  2.2,  3.7,  and  4.7  pm  regions.  Figure  4  shows  typical  ARES  75-channel  spectra 
obtained  over  clear  and  cirrus  cloudy  areas.  Respective  spectral  intervals  for  the  absorption  bands  are 
noted  in  the  diagram.  Large  measurement  noises  exist  for  Chs.  1-25  (1. 9-2.7  pm),  where  the  solar 
reflection  contributions  are  dominant,  and  the  cloud  reflection  is  stronger  than  the  surface  reflection.  For 
Chs.  26-75,  the  aforementioned  four  absorption  bands  and  two  window  bands  are  clearly  displayed  as 
local  minima  and  maxima  within  both  spectra,  respectively.  Over  the  3.7  pm  and  4.7  pm  window  bands 
and  the  6.3  pm  water  vapor  absorption  band-wing  interval,  there  are  significant  differences  between  the 
clear  and  cloudy  radiances.  These  differences  can  be  used  to  detect  the  presence  of  cirrus  clouds  and  to 
retrieve  their  optical  and  microphysical  parameters. 

3.2  Selection  of  Retrieval  Channels 

To  select  a  number  of  appropriate  channels  for  the  cirrus  cloud  retrieval,  we  have  examined  all 
the  ARES  channel  weighting  function  profiles  and  images.  Typical  weighting  function  profiles  were 
generated  using  the  vertical  transmittance  functions  obtained  from  the  MODTRAN  calculations 
employing  the  climatological  midlatitude  summer  temperature  and  humidity  soundings.  These  weighting 
function  profiles  show  features  that  correspond  to  the  distribution  of  absorption  bands.  For  the  window 
channels  (3.75  pm,  and  4.7  pm),  the  weighting  function  peaks  at  or  near  the  surface,  so  that  the  clear 
radiance  is  dominated  by  surface  emission.  The  inhomogeneity  of  surface  types  and  emissivity  can 
cause  clear  radiance  to  spread,  leading  to  errors  in  cloud  retrieval  results.  For  the  2.7  pm  H2O,  4.3  pm 
CO2.  and  6.3  pm  H2O  band-center  channels,  the  weighting  function  peaks  in  the  upper  atmosphere.  The 
radiance  tends  to  be  dominated  by  gaseous  emission  from  the  upper  atmosphere,  suppressing  the  cloud 
signals.  For  the  6.3  pm  band-wing  channels,  the  weighting  function  peaks  in  the  middle  and  lower 
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troposphere.  For  the  purpose  of  the  detection  and  retrieval  of  cirrus  clouds,  it  is  advantageous  to  choose 
channels  with  the  weighting  functions  peaking  in  the  lower  troposphere.  This  is  because  the  data  spread 
in  clear  radiance  is  reduced  due  to  the  relatively  homogeneous  and  moderate  water  vapor  emission  in  the 
band-wing  spectral  interval,  and  the  upwelling  emission  of  cirrus  clouds  can  be  unambiguously  detected 
by  satellite  radiometers  without  attenuation  by  water  absorption/emission  in  the  upper  atmosphere. 

Moreover,  we  have  used  radiance  images  to  facilitate  the  selection  of  channels  by  acquiring  the 
ARES  data  collected  on  September  16,  1995  (No.  16-9-95-1  A).  Using  a  quick-look  software  ADAPT 
(ARES  Data  Analysis  and  Processing  Toolkit  developed  by  the  SCITEC  Inc.  at  Princeton,  N.J),  we  can 
display  graphically  the  large  volume  of  the  ARES  dataset  on  a  SGI  system  in  imagery  forms.  We  have 
manually  e.xamined  the  image  of  each  flight  track  for  each  ARES  channel,  and  have  selected  a  scene 
which  displays  a  distinct  contrast  between  high-  and  low-  radiance  areas  using  the  4.747  pm  channel 
imagery.  Based  on  both  ground-based  8.6  mm  radar  and  airborne  Ka-band  radar  and  0.53  pm  lidar 
observations  on  board  HARP,  only  a  single  layer  of  cirrus  cloud  was  present  during  the  flight  period. 
Thus,  the  high-  and  low-radiance  areas  correspond  to  clear  and  cloudy  regions,  respectively.  Figure  5 
shows  the  radiance  images  for  Chs.  41(3.681  pm).  42(3.755  pm),  55(4.747  pm),  60(5.116  pm),  and 
63(5.332  pm)  over  the  selected  scene.  The  first  three  channels  are  within  the  3.75  pm  and  4.7  pm 
window  bands,  while  the  last  two  channels  are  within  the  6.3  pm  water  vapor  band-wing  interval.  For 
each  image,  the  horizontal  and  vertical  scales  are  -500  m  and  2.5  km,  respectively. 

From  these  images,  it  is  noted  that  the  upper  1/3  of  the  domain  is  clear,  while  the  rest  of  the 
domain  is  cloudy.  Over  the  clear  area,  the  images  for  the  three  window  channels  show  significant 
variations  of  the  emitted/reflected  radiances  due  to  the  inhomogeneity  of  the  surface  temperature  and 
emissivity.  However,  for  the  two  band-wing  channels,  radiances  over  the  clear  area  appear  to  be 
relatively  uniform,  as  compared  to  the  clear  radiances  for  the  window  channels,  because  the  significant 
water  vapor  absorption/emission  effect  in  the  lower  troposphere  covers  the  horizontally  inhomogeneous 
surface  emission.  The  narrow  distribution  of  the  clear  radiances  is  a  desirable  feature  for  application  of 
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the  remote  sensing  algorithm.  A  smaller  spread  in  clear  radiances  would  lead  to  less  errors  in  prescribing 
the  mean  clear  radiances,  thus  increasing  the  retrieval  accuracy  (Rao  et  al.  1995). 

Over  the  cloudy  area,  it  is  evident  that  the  image  for  the  4.7  pm  channel,  as  well  as  the  images  for 
the  two  band-wing  channels,  contain  distinct  cloud  signatures.  For  the  two  3.7  pm  window  channels, 
there  is  no  distinguishable  cloud  signature,  because  the  compensating  effects  of  the  strong  reflection  and 
weak  emission  by  cloud  particles  cause  the  resulting  total  cloud  radiance  to  be  about  the  same  as  the  total 
clear  radiance.  Radiance  images  for  channels  at  the  absorption  band  center  (not  shown  here)  exhibit  no 
cloud  signature.  Based  on  the  preceding  discussions  on  the  properties  of  the  computed  weighting 
function  profiles  and  the  images  shown  in  Figure  5,  we  select  Ch.  60.  61  (5.189  pm),  62  (5.261  pmj  and 
63  for  the  development  of  a  cirrus  cloud  retrieval  algorithm.  The  weighting  function  for  these  channels 
peaks  in  the  lower  troposphere,  and  images  for  these  channels  contain  definitive  cloud  signatures  and 
display  narrow  data  spread  of  clear  radiances. 

3.3  Retrievals  using  ARES  Channel  Data 

In  order  to  apply  the  aforementioned  algorithm  to  the  inference  of  the  cloud  temperature  and 
emissivity,  we  must  first  determine  clear  or  cloudy  pixels  within  any  selected  scene.  Thus,  we  have 
developed  a  clear/cloud  detection  scheme  to  identify  clear  and  cloudy  data  points.  It  is  based  on  the 
physical  properties  of  the  ARES  equivalent  brightness  temperatures  for  selected  channels  for  clear  and 
cloudy  pixels.  In  this  scheme,  clear  pixels  are  separated  from  the  cloudy  pixels  by  applying  the  threshold 
tests  to  the  equivalent  brightness  temperature  data  for  each  selected  channels.  If  the  equivalent  brightness 
temperature  for  a  pixel  is  greater  than  the  prescribed  threshold  value  for  each  selected  channels,  then  it  is 
categorized  as  clear.  Otherwise,  it  is  flagged  as  cloudy.  The  resulting  number  of  cloudy  pixels  may  be 
overestimated  in  this  scheme.  However,  by  doing  so,  we  can  ascertain  that  the  detected  clear  pixels  have 
the  maximum  probability  to  be  truly  clear,  and  that  the  maximum  amount  of  cloudy  pixels  are  available 


for  the  retrieval  to  be  carried  out. 
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Appropriate  threshold  values  for  the  selected  channels  are  determined  from  statistical  analyses  on 
the  radiance  data  for  the  particular  case  of  the  mission  flight  on  September  16,  1995.  Figures  6(a)-(d) 
show  the  histograms  of  the  brightness  temperatures  for  Chs  60-63  within  the  domain  of  Figure  5.  Each 
histogram  displays  a  distinct  bi-modal  shape  due  to  the  presence  of  contrasting  clear  and  cloudy  areas 
within  the  selected  scene.  As  discussed  previously,  the  high-  and  low-radiance  modes  represent  the  mean 
radiances  for  the  clear  and  cloudy  conditions,  respectively.  From  these  analyses,  the  threshold  value  for 
each  channel  can  be  easily  determined  as  the  radiance  that  is  associated  with  the  minimum  frequency 
between  the  two  modes.  Once  the  clear  pixies  are  defined,  the  mean  clear  radiance  for  each  channel  is 
determined  based  on  the  average  of  radiances  of  all  detected  clear  pixels  within  the  selected  scene  (-10 
km  stripi. 

Next  we  investigate  the  characteristics  of  the  measured  ARES  channel  data.  The  ARES  thermal 
IR  channel  data  exhibit  characteristics  of  linearity  in  the  two-channel  radiance  correlation  as  implied  in 
Eq.  (3).  Figure  7  shows  the  scatter  diagram  for  the  three  pairs  of  ARES  channels:  61-60,  62-60,  and  63- 
60.  In  each  frame,  the  curve  which  is  nearly  linear  denotes  the  correlation  of  Planck  functions  for  the  two 
channels,  and  the  dots  are  radiance  data  points  chosen  from  the  scene  of  Fig.  5  covering  the  area  of  the 
strong  contrast  between  clear  and  cloudy  regions.  The  manner  in  which  radiance  data  points  distribute  in 
clusters  indicates  a  statistically  significant  linear  correlation  between  each  pair  of  channel  radiances. 
There  are  two  clusters  in  each  frame.  The  upper  cluster  corresponds  to  clear  data  points,  while  the  lower 
elongated  clusters  are  associated  with  cloudy  radiances. 

Figure  8  demonstrates  the  comparability  of  the  measured  and  computed  3.7  pm  channel  radiances 
for  the  retrieval  of  mean  effective  sizes.  In  the  upper  panel,  the  computed  3.7  pm  solar  radiances  are 
plotted  as  functions  of  optical  depth  and  six  mean  sizes.  For  T  >  2,  these  radiances  are  nearly  independent 
of  the  optical  depth,  but  are  strong  functions  of  the  ice  crystal  size  distribution.  In  the  middle  panel,  the 
computed  3.7  pm  thermal  IR  radiances  for  different  sizes  are  plotted  against  the  optical  depth.  In  this 
case,  the  cloud  temperature  is  250  K  and  the  clear  brightness  temperature  is  291  K.  The  3.7  pm  thermal 
IR  radiances  decrease  with  increasing  optical  depth  and  mean  effective  size.  Finally,  in  the  lower  panel, 
the  computed  3.7  pm  total  radiances  are  plotted  against  the  optical  depth  for  the  six  ice  crystal  sizes.  It  is 
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noted  that  these  radiances  are  sensitive  to  both  the  optical  depth  and  ice  crystal  size.  We  display  the 
ARES  measured  data  collected  during  the  Track  3  of  the  September  16  flight  on  the  lower  diagram  and 
show  that  they  lie  within  the  bulk  range  of  the  computed  curves,  indicating  that  it  is  feasible  to  use  the 
ARES  measured  data  to  retrieve  ice  mean  effective  sizes. 

Applying  the  retrieval  scheme  to  the  detected  cloudy-pixel  radiance  data  of  the  selected  ARES 
channels,  we  can  determine  the  cirrus  cloud  temperatures  and  emissivities  (optical  depth)  for  the  four 
flight  tracks  of  September  16,  1995.  By  using  an  efficient  numerical  scheme,  the  total  run  time  on  a 
powerful  state-of-the-art  SGI  is  less  than  10  min  for  each  track.  For  the  purpose  of  displaying  data,  we 
obtain  the  array  averaged  channel  equivalent  brightness  temperatures  and  the  retrieved  parameters  for 
each  scan  within  the  flight  track.  Variation  of  the  computed  array-averaged  Ch.  60  equivalent  brightness 
temperatures,  cloud  temperatures,  optical  depths,  and  mean  effective  sizes  for  Track  1  is  shown  in  Fig.  9 
based  on  the  application  of  the  retrieval  scheme  to  10"*  scans  of  the  45-pixel  array  data.  There  are  three 
periods  which  have  been  identified  as  clear  (-16-20  km,  25-28  km  and  30-32  km).  The  retrieved  cloud 
temperatures,  optical  depths,  and  mean  effective  ice  crystal  sizes  vary  between  230  and  260  K,  between  0 
and  2,  and  between  25  and  75  [im,  respectively.  Further  examination  of  these  retrieval  results  reveals  that 
the  cirrus  cloud  optical  depths  were  generally  less  than  1.0,  while  the  cirrus  cloud  temperatures  were 
mostly  higher  than  -30°C.  The  cirrus  clouds  were  found  to  be  composed  of  relatively  small  ice  crystals. 
It  is  possible  that  these  cirrus  clouds  were  undergoing  through  a  formation  stage,  during  which  ice  crystal 
sizes  are  usually  small  and  clouds  are  optically  thin.  The  reliability  of  the  aforementioned  retrieval  results 
has  been  confirmed  by  the  following  validation  efforts. 

3.4  Validation  of  Retrieval  Results 

To  validate  the  retrieved  cirras  cloud  temperature/height,  we  use  the  time  series  of  cloud 
boundary  (cloud  top  and  base)  altitudes  derived  from  the  8.6-mm  radar  measurements,  together  with  the 
upper  tropospheric  temperature  measurement  by  the  Rosemount  total  temperature  sensor  on  board  the 
HARP.  Figure  10(a)  shows  the  atmospheric  temperature  at  various  levels  between  5  and  11  km.  These 
temperature  values  were  obtained  between  2029  and  2341  UTC.  They  vary  between  -5°  and  -55°C 
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within  the  altitude  range.  This  temperature  profile  is  used  in  the  validation  effort  because  there  was  no 
surface  sounding  launched  around  the  flight  time  of  the  WB-57F.  We  assume  that  the  upper  atmospheric 
thermal  structure  remained  the  same  during  the  HARP  flight  period,  so  that  the  temperature  profile  shown 
in  Fig.  10(a)  is  representative  throughout  this  period.  Figure  10(b)  shows  the  temporal  variation  of  the 
cirrus  cloud  boundary  altitudes  derived  from  the  8.6-mm  radar  measurements  between  2000  and  2400 
UTC.  A  single  layer  of  cirrus  cloud  was  present  during  this  time  period.  Between  2(XX)  and  2100  UTC, 
the  observed  cloud  top  and  base  were  located  at  8  and  6  km,  respectively.  After  2100  UTC,  the  cloud-top 
altitude  fluctuates  between  8  and  10  km,  while  the  cloud  base  altitudes  largely  remains  at  6  km.  except  for 
the  period  between  2130  and  2230  UTC.  There  are  two  clear/thin-cirrus  episodes  around  2100  and  2140- 
2200  UTC  and  one  broken-cirrus  period  around  2210-2220  UTC.  The  retrieved  cloud  temperatures  were 
converted  into  cloud  altitudes  using  the  temperature  sounding  shown  in  Fig.  10(a).  For  Tracks  1  and  3, 
when  WB-57F  flew  over  the  radar  site,  the  area  happened  to  be  clear,  although  the  radar  recorded  very 
thin  cirrus.  Thus,  comparison  of  the  retrieved  cloud  altitudes  with  the  radar  data  cannot  be  accomplished. 
For  Tracks  5  and  8,  the  retrieved  cloud  temperatures  over  the  radar  site  are  -33°  and  -40  °C,  which 
correspond  to  8.2  and  9.2  km,  respectively  (denoted  by  square  boxes  in  Fig.  10(b)).  Both  altitudes  are 
within  the  radar-derived  cloud  boundaries  but  near  the  cloud  top  at  the  particular  WB-57F  overpass  time. 

To  verify  the  retrieved  cirrus  cloud  optical  depths  and  ice  crystal  mean  effective  sizes,  we  use  the 
in-situ  ice  crystal  size  distributions  obtained  by  the  optical  probes  on  board  the  HARP  to  compute  the 
single-scattering  properties  of  the  ice  crystals  and  to  evaluate  the  in-situ  mean  effective  size.  A 
representative  ice  crystal  size  distribution  has  been  constructed  for  each  case  by  using  the  2D-C  probe 
measurements  for  ice  crystals  larger  than  75  jim,  and  an  extrapolation  technique  for  ice  crystals  smaller 
than  75  |J.m  in  a  manner  described  in  the  following. 

Among  the  cloud  microphysic  instruments  on  board  HARP,  the  FSSP  sizes  the  ice  crystals  from 
1  to  96  pm;  the  1 D-C  probe  samples  ice  crj'stals  from  20  to  300  pm;  and  the  2D-C  probe  measures  the 
maximum  dimensions  of  the  projection  of  ice  crystals  from  25  to  800  pm.  The  size  distributions  derived 
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from  the  ID-C  and  2D-C  probes  closely  agree  with  each  other  for  the  size  range  of  25  to  300  ^m.  though 
the  differences  in  the  ice  crystal  number  concentrations  for  some  size  bins  can  be  as  large  as  one  order  of 
magnitude.  We  choose  to  use  the  2D-C  derived  ice  crystal  size  distribution  because  this  probe  covers  the 
typical  ice  crystal  size  range  and  because  the  FSSP  derived  size  distributions  appear  to  be  unreliable  in 
view  of  the  unrealistically  high  ice  crystal  number  concentrations  for  sizes  less  than  20  pm.  For  these 
reasons,  we  use  the  Gamma  function  to  extrapolate  the  2D-C  size  distribution  down  to  1  pm  to  obtain  the 
ice  crystal  size  distribution  for  ice  crystal  sizes  less  than  75  pirn  (Tsay  et  al.  1996).  Composite  ice  crystal 
size  distributions  combining  the  2D-C  probe  measurements  and  the  extrapolated  portions  are  plotted  in 
Fig.l  1  (a)  and  (b).  The  upper  bound  for  both  the  ice  ciy  stal  size  distributions  is  less  600  pm.  Both  ice 
crystal  size  distributions  are  similar  to  those  occurred  in  typical  cirrostratus  clouds  (Heymsfield  1975). 
We  find  that  incorporating  the  extrapolated  portion  from  20  pm  down  to  1  pm  has  negligible  effects  on 
the  extinction  coefficient.  Other  single-scattering  parameters  are  not  sensitive  to  the  inclusion  of  ice 
crystals  smaller  than  75  pm  for  the  size  distributions  used  in  the  present  smdy.  Visible  extinction 
coefficients  of  0.54  and  0.57  km''  were  obtained  for  cases  1  and  2.  respectively. 

To  determine  the  optical  depth  based  on  in-situ  measurements,  the  average  cirrus  cloud  thickness 
and  the  associated  extinction  coefficients  must  be  known.  Based  on  the  cloud  boundary  altitudes  derived 
from  the  8.6-mm  radar  measurements,  the  averaged  cirrus  cloud  thicknesses  during  the  two  periods  were 
about  1 .75  and  2  km  for  cases  1  and  2,  respectively.  The  ascending  vertical  profile  of  the  2D-C  median 
volume  ice  crystal  diameter  indicates  that  ice  crystal  sizes  vary  by  no  more  than  20  pm  between  6.5  and  8 
km  (Morrison  et  al.  1997).  For  this  reason,  we  use  the  computed  extinction  coefficients  based  on  single- 
level  in-situ  size  distributions  to  represent  the  vertically  averaged  cirrus  cloud  extinction  coefficients. 
The  computed  optical  depths  based  on  the  preceding  approximation  are  0.94  and  1.14  for  cases  1  and  2, 
respectively.  The  mean  retrieved  optical  depths  are  0.89  and  1.42  for  Tracks  1  and  3  of  the  ARES  flight, 
respectively.  The  retrieved  mean  optical  depth  for  Track  1  agrees  relatively  well  with  that  from  the  2D-C 
probe.  However,  for  Track  3,  the  difference  between  the  two  optical  depths  is  on  the  order  of  0.3. 
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Finally,  we  utilize  a  cirrus  cloud  model  composed  of  50%  bullet  rosettes,  30%  hollow  columns, 
and  20%  plates  to  compute  the  in-situ  mean  effective  ice  crystal  sizes.  An  empirical  relationship  between 
the  maximum  dimension  (L)  and  width  (W)  is  employed  for  each  shape;  for  plates,  W  =  2.02,  L  for 
hollow  columns,  W  =  1 1.3,  L  for  L  >  200  pm,  W  =  8.479  +  1.002  L  -  0.00234  L  ^  for  L  <  200  pm; 
and  for  bullet  rosettes,  W  =  2.3103  The  in-situ  mean  effective  sizes  are  based  on  the  average  of  the 
mean  effective  size  for  each  shape  weighted  by  the  percentage  distribution.  They  are  57.2  pm  for  case  1 
and  55.2  pm  for  case  2.  The  retrieved  mean  effective  sizes  are  64.2  pm  for  Track  1  and  53.3  pm  for 
Track  3,  which  compare  reasonably  well  with  the  in-situ  sizes  with  differences  less  than  10  pm.  There  are 
two  reasons  for  differences  between  retrieved  and  observation-derived  optical  depths  and  mean  effective 
sizes.  First,  estimate  of  the  optical  depth  from  the  2D-C  probe  is  based  on  the  composite  ice  crv'stal  size 
distribution  at  a  single  flight  level  within  the  cirrus  cloud.  Second,  the  variation  range  of  the  retrieved 
optical  depth  for  Track  3  is  much  larger  than  that  for  Track  1. 

4.  Summary 

We  have  presented  a  remote  sensing  algorithm  for  the  retrieval  of  cirrus  cloud  temperature 
(height),  optical  depth,  and  ice  crystal  size  using,  for  the  first  time,  the  spectral  channels  in  the  5.1 -5.3pm 
range  in  combination  with  the  3.7  pm  channel.  This  scheme  has  been  applied  to  the  ARES  data  collected 
on  September  16,1995,  over  the  western  Boston  area.  We  have  discussed  the  ARES  channel  spectral 
properties  and  the  rationale  for  the  selection  of  retrieval  channels.  Characteristics  of  the  measured  ARES 
channel  data  have  been  examined  to  demonstrate  their  suitability  for  retrieval.  For  retrieval  purposes,  a 
simple  clear/cloud  threshold  detection  scheme  has  been  developed  and  applied  to  identify  clear  and 
cloudy  data  points. 

Validation  of  the  retrieved  cirrus  cloud  temperature/height  has  been  carried  out  by  using  the 
collocated  and  coincident  time  series  of  the  cloud  boundary  (cloud  top  and  base)  altitudes  derived  from 
the  ground-based  8.6  mm  radar  measurements.  We  have  also  compared  the  retrieved  cirrus  cloud  optical 
depths  and  ice  crystal  mean  effective  sizes  with  those  parameters  determined  based  on  the  composite  ice 


15 


crystal  size  distributions  from  in-situ  measurements.  We  show  that  the  retrieved  cloud  heights  compare 
reasonably  well  with  the  values  determined  from  the  ground-based  8.6-mm  radar,  and  that  the  retrieved 
optical  depth  and  mean  effective  size  match  the  values  determined  from  the  in-situ  microphysical 
measurements.  Differences  between  the  retrieved  and  observation-derived  optical  depths  and  mean 
effective  sizes  could  be  due  to  the  limitation  in  data  collection  by  aircraft  and  the  inherent  cloud 
horizontal  and  vertical  inhomogeneity. 
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FIGURE  CAPTIONS 

The  k-value  as  a  function  of  De  for  the  retrieval  channels.  Eleven  size  distributions  have  been 
used  for  single-scattering  computations.  The  single-scattering  properties  are  used  in  an 
adding/doubling  code  to  obtain  the  cloud  emissivity. 

A  schematic  flow  chart  for  the  cirrus  retrieval  algorithm  using  ARES  channel  data. 

Upwelling  radiances  for  solar  reflection  and  thermal  emission  in  the  2-6.5  pm  region.  F 
denotes  the  solar  irradiances,  p  is  the  cosine  of  solar  zenith  angle,  a  is  the  surface  albedo,  and 
Bx  is  the  Planck  function. 

Typical  ARES  75-channel  spectra  for  clear  and  cloudy  areas. 

Radiance  images  for  ARES  Chs.  41,  42,  55,  60,  and  63  over  a  selected  scene  having  a  strong 
cloud/clear  contrast  within  the  Track  1  of  September  16  flight. 

Histograms  of  the  brightness  temperatures  for  the  ARES  Chs.  60-63  within  the  domain  of  Fig. 
7. 

Examples  of  the  scatter  diagram  for  three  pairs  of  ARES  channels:  60-61,  60-62,  and  60-63. 
Computed  3.7  pm  solar  radiances  (a),  IR  radiances  (b),  and  total  radiances  (c)  as  functions  of 
optical  depth  and  ice  crystal  size. 

Array  averaged  channel  equivalent  brightness  temperatures,  retrieved  cloud  temperature,  cloud 
optical  depth,  and  effective  size  as  functions  of  flight  distance  for  Track  I,  9/16/95. 

(a)  Partial  sounding  of  the  middle  and  upper  troposphere,  obtained  from  the  HARP  flight  above 
Hanscom  AFB  (squares  correspond  to  the  retrieved  cloud  temperature),  (b)  Comparison  of  the 
cloud  top  heights  obtained  from  radar/lidar  measurements  and  the  retrieved  cloud  temperatures. 
Retrieved  temperatures  are  displayed  for  Tracks  1,  3,  5,  and  8  (in  chronological  order). 
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Fig.  11.  The  ice  crystal  size  distributions  based  on  the  HARP  data  (9/16/95)  with  the  lower  limit  of  75 
urn  and  extrapolations  to  20  pm  and  1  pm  based  on  the  Gamma  function  and  an  optimization 
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